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This  dissertation  presents  a  quantitative  study  of  the  physical 
mechanisms  underlying  the  anomolously  large  recombination  current 
experimentally  observed  in  heavily  doped  regions  of  silicon  pn-junction 
solar  cells  and  bipolar  transistors.  The  study  includes  a  comparison 
of  theoretical  predictions  with  a  variety  of  experimental  observations 
in  heavily  doped  silicon  and  silicon  devices. 

A  major  conclusion  is  that  the  simplest  physical  model  that 
adequately  describes  the  heavily  doped  regions  must  include  Fermi- 
Dirac  statistics,  a  phenomenological  excess  intrinsic  carrier  density 
(or  deficit  impurity  concentration),  Auger  recombination  in  the  bulk, 
and  recombination  at  the  surface.  These  mechanisms  are  incorporated  in 
a  first-order  model  useful  in  the  design  of  silicon  pn-junction  solar 
cells.  The  accuracy  of  the  first-order  model  is  supported  by  comparing 
its  results  with  the  results  of  more  detailed  models  and  of  a  numerical 


vn 


analysis  of  the  problem.  Experimental  data  are  presented  that  are 
consistent  with  the  predictions  of  the  first-order  model  and  of  the 
numerical  solution. 
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CHAPTER  I 
INTRODUCTION 

Heavily  doped  regions  of  silicon  pn-junction  devices  show 
experimentally  a  larger  recombination  current  than  is  predicted  by  classi- 
cal pn-junction  theory  [1,2].  This  large  recombination  current  in  heavily 
doped  regions  is  responsible  for  the  low  values  of  the  open-circuit 
voltage  VQC  of  silicon  solar  cells  [3]  and  the  common-emitter  current 
gain  hpE  in  silicon  bipolar  transistors  [4,5]. 

The  main  objective  of  this  dissertation  is  to  study  quantitatively 
large  recombination  currents  in  heavily  doped  regions  and  to  describe 
the  fundamental  limitations  imposed  by  heavy  doping  on  the  performance 
of  silicon  pn-junction  solar  cells  and  bipolar  transistors.  The 
achievement  of  this  objective  involves: 

A.  Identifying  physical  mechanisms  associated  with 
heavy  doping. 

B.  Determining  the  dominant  heavy-doping  mechanisms 
that  affect  the  device  performance. 

C.  Including  the  dominant  heavy-doping  mechanisms  in 
analytic  and  computer-aided  analysis  of  heavily 
doped  regions  in  silicon  devices. 

D.  Establishing  accurate  and  simple  first-order 
engineering  design  models  based  on  device  physics. 


From  charge-control  theory  [6-10],  the  recombination  current 
in  any  heavily  doped  region,  for  example,  the  minority-carrier  emitter 
current  J^  can  be  expressed  by 

J  -Si 

E   TE  (1.1) 

where  QE  is  the  minority-carrier  charge  storage  and  te  is  the  effective 
lifetime  in  the  emitter.  From  classical  pn-junction  theory  [11],  QE 
increases  with  decreasing  bandgap  and  decreases  with  increasing  majority 
carrier  density.  The  effective  lifetime  ~E  may  be  the  average  bulk 
lifetime  of  the  minority  carriers  or  their  transit  time,  defined  as  the 
average  time  needed  by  a  minority  carrier  to  cross  the  semiconductor 
region  without  recombining  in  the  bulk,  or  a  combination  of  the  two. 
From  (i.l)  we  can  see  that  an  excessive  emitter  current  is  due  to  (a)  a 
large  minority-carrier  charge,  (b)  a  short  effective  lifetime,  and 
(c)  a  combination  of  (a)  and  (b). 

Various  physical  mechanisms  associated  with  heavy  doping  in 
semiconductors  can  give  rise  to  (a),  (b),  or  (c).  These  heavy-doping 
mechanisms  are  outlined,  in  Table  I,  and  contrasted  to  corresponding 
mechanisms  in  lightly  doped  semiconductors.  We  will  describe  briefly 
some  of  the  fundamental  heavy-doping  mechanisms  outlined  in  Table  I. 

One  of  the  fundamental  heavy-doping  mechanisms  that  yields  a 
large  minority-carrier  charge,  and  consequently  a  large  current  by 
(1.1)?  is  energy-bandgap  narrowing.  From  a  theoretical  standpoint, 
various  mechanisms  exist  that  lead  to  bandgap  narrowing  [12-19].  More- 
over, electrical  measurements  in  heavily  doped  single-crystal  silicon 


TABLE  I 

COMPARISON  BETWEEN  A  LOW  IMPURITY  CONCENTRATION 

SEMICONDUCTOR  (LESS  THAN  a  10"16/CM3)  AND  A  HIGH 

IMPURITY  CONCENTRATION  SEMICONDUCTOR. (MORE  THAN  *  10^/CM3] 


Low  Impurity  Concentration 


High  Impurity  Concentration 


Effective  bandgap =  Intrinsic  band-  Effective  bandgap < Intrinsic  band- 
gap  gap 


(AEQ  =  0) 


(aeg^o; 


Maxwell -Boltzmann  statistics 
for  the  majority  and  minority 
carriers 


Maxwell -Boltzmann  statistics  for 
minority  carriers;  Fermi-Dirac 
statistics  for  majority  carriers 


Simple  Law  of  Mass  Action  valid    Generalized  Law  of  Mass  Action 


P0N0  =  n2(T) 


P0N0  =  nj  (T,x) 


All  impurity  atoms  ionized  (at 
room  temperature) 


Seme  impurity  atoms  deionized  (at 
room  temperature) 


V*\,   Up,  Dp  and  D^  constant  (in-    u^,  yp,  Dp  and  Dm  dependent  on 
dependent  of  impurity  concentra-    impurity  concentration 
tion) 


Einstein's  relation  valid 
D  _kT 

u  q 


Einstein's  relation  not  valid; 
modified  relation: 


1/2 


(n! 


D  x  JcT  ,. 

y  "  q      F-i/2  ^n' 


Built-in  electric  field  in  QNR 
is  same  for  majority  and  minor- 
ity carriers 


Effective  drift  field  in  QNR  is  not 
same  for  majority  and  minority 
carriers 


maj   mm 


E   •  r  E  • 
Lmaj  r  cmin 


TABLE  I--Continued 


Low  Impurity  Concentration 


High  Impurity  Concentration 


Recombination  of  excess  minor- 
ity carriers  via  defects, 
(Shockley-Read-Hall) 


Recombination  of  excess  minority 
carriers  via  (a)  fundamental  band- 
to-band  Auger  process,  (b)  position 
dependent  recombination  centers  or 
defects  (Shockley-Read-Hall) 


materials  and  devices  indicate  that  the  forbidden  energy  gap  may  be 
smaller  than  the  intrinsic  energy  gap  of  pure  silicon.  These  experi- 
mental measurements  are: 

A.  Measurements  of  the  temperature  dependence  of  the 
current  gain  in  silicon  bipolar  transistors  [4,5] 

B.  Measurements  of  the  collector  current  and  its 
temperature  dependence  in  silicon  bipolar  tran- 
sistors [20] 

C.  Measurements  of  the  temperature  dependence  of  the 
photoresponse  of  silicon  pn-junctions  and  tran- 
sistors [21] 

D.  Measurements  of  the  temperature  dependence  of  the 
emitter  current  in  silicon  solar  cells  (in  the  dark 
condition)  and  bipolar  transistors  [22] 

E.  Measurements  of  the  emitter  current  in  different 
pn-junction  structures  [23] 

Rather  than  providing  a  direct  measurement  of  the  bandgap,  the 

above  measurements  yield  an  excess  effective  intrinsic  carrier  density 

2 
nie,  defined  as  the  thermal  equilibrium  pn-product  in  heavily  doped 

silicon  [24].  The  measured  values  of  n^  are  considerably  larger  than 

n.   for  pure  silicon. 

Another  mechanism  that  can  provide  an  excess  intrinsic  carrier 

density,  and  a  large  minority-carrier  charge  in  the  heavily  doped 

region,  is  deionization  of  impurities  at  high  doping  concentrations 

[25].  For  low  concentrations  of  shallow  impurities  in  silicon  at  room 

temperature,  all  impurity  atoms  are  essentially  ionized.  This  is  so 


because  the  Fermi  level  is  well  below  the  impurity  levels.  As  the 
impurity  concentration  increases  (above  1016  cm"3)  some  deionization 
of  impurities  occur.  The  degree  of  ionization  of  impurity  atoms  depends 
on  the  density  and  position  of  localized  energy  levels  in  the  bandgap 
and  on  the  position  of  the  Fermi  level  relative  to  these  levels. 

Another  fundamental  heavy-doping  mechanism  that  must  be  con- 
sidered is  the  degeneracy  of  the  majority  carriers  in  heavily  doped 
semiconductors.  Majority-carrier  concentrations  in  heavily  doped 
semiconductors  are  high  enough  that  the  carriers  can  interact  with 
each  other  via  the  Pauli  Exclusion  Principle.  Thus  Maxwell -Boltzmann 
statistics  are  no  longer  applicable  and  Fermi -Dirac  statistics  must  be 
used  instead. 

We  consider  now  two  other  fundamental  heavy-doping  mechanisms 
that  can  increase  J£  in  (1.1)  by  decreasing  te.  The  first  mechanism 
is  Auger  band-to-band  recombination  [26].  It  is  fundamental  in  the 
sense  that  it  sets  an  upper  bound  on  the  minority-carrier  lifetimes 
in  these  regions.  The  second  mechanism  is  Schockley-Read-Hall 
recombination  [27-28]  via  position-dependent  defect  densities. 

We  have  discussed  the  fundamental  heavy-doping  mechanisms 
that  can  limit  the  performance  of  heavily  doped  silicon  pn-junction 
devices  such  as  solar  cells  and  bipolar  transistors.  We  outline  now 
our  approach  in  assessing  the  relative  importance  of  these  heavy-doping 
mechanisms  to  tne  performance  of  silicon  devices. 

Our  approach  in  assessing  the  importance  of  heavy-doping 
effects  is  to  include  each  of  these  effects  separately  in  a  simple  yet 


accurate  model  for  the  heavily  doped  region.  He  first  include  Auger 
recombination  alone,  neglecting  other  heavy-doping  effects,  in  a 
rigorous  analytic  model  for  the  emitter  region  of  silicon  solar  cells 
and  bipolar  transistors.  The  model  includes  also  a  finite  surface 
recombination  velocity  at  the  emitter  surface.  We  find,  in  Chapter  II, 
that  values  of  VqC  in  silicon  solar  cells,  predicted  by  this  model, 
are  considerably  higher  than  those  experimentally  observed.  So 
Auger  recombination  is  not  the  dominant  mechanism  that  yields  the  large 
emitter  current  in  bipolar  devices. 

Next,  we  consider  recombination  via  defects,  in  Chapter  III. 
Using  a  previous  detailed  study  of  defect  recombination  [9],  we 
argue  that  if  recombination  via  defects  is  the  main  mechanism  underlying 
the  large  emitter  current,  then  surface  recombination  should  not 
influence  the  device  performance.  This  implication  is  inconsistent  with 
recent  experimental  data  [29]  showing  the  sensitivity  of  Vqq  in  typical 
silicon  solar  cells  to  surface  treatment.  These  experimental  data 
emphasize  the  importance  of  surface  recombination  and  an  excess  in- 
trinsic  carrier  density  in  heavily  doped  silicon  (i.e.,  r\.     >   n^). 

The  importance  of  deionization  in  heavily  doped  silicon  is 
assessed  in  Chapter  IV.  We  investigate  whether  deionization  can  be  the 
main  mechanism  underlying  the  experimentally  observed  effective 
intrinsic  carrier  density  in  heavily  doped  silicon.  By  using  a 
variety  of  experimental  data,  we  infer  that  deionization  in  heavily 
doped  silicon  at  room  temperature  is  not  an  important  mechanism,  at  least 
not  in  accordance  with  the  recent  models  of  Pooovic  and  Heasell . 


Having  established  the  important  fundamental  heavy-doping 
mechanisms,  we  proceed  to  develop  engineering  design  models  for 
heavily  doped  regions. 

Energy-bandgap  narrowing,  Fermi-Dirac  statistics,  and  a  finite 
surface  recombination  velocity  are  included  in  a  rigorous  analytic 
model  of  the  heavily  doped  emitter  of  silicon  devices  in  Chapter  V. 
The  emitter  is  assumed  to  be  transparent  to  the  injected  minority 
carriers,  that  is,  most  of  the  injected  minority  carriers  can  cross 
the  quasi-neutral  emitter  region  and  recombine  mainly  at  the  surface 
rather  than  in  the  bulk.  The  transparency  assumption  can  be  tested 
for  self-consistency  for  any  given  device.  He  show,  in  Chapter  V. 
that  the  transparent  emitter  model  can  explain  the  experimentally 
observed  values  of  VQC  in  silicon  solar  cells  if  the  surface  recombina- 
tion velocity  is  high  and  energy-bandgap  narrowing  and  Fermi-Dirac 
statistics  are  included.  In  agreement  with  other  experimental  data, 
we  also  show  that  reducing  the  front  surface  recombination  velocity 
increases  Vqc  in  silicon  solar  cells. 

In  Chapter  VI,  we  develop  a  simple  first-order  model  for  the 
emitter  that  includes  energy-bandgap  narrowing,  Fermi-Dirac 
statistics,  Auger  recombination,  and  a  finite  surface  recombination 
velocity.  The  first-order  model  is  useful  as  a  design  tool  particularly 
for  silicon  solar  cells.  The  accuracy  of  the  model  is  checked  by  a 
rigorous  computer-aided  numerical  analysis  of  the  problem.  The  model 
is  found  to  be  accurate  if  the  surface  impurity  concentration  is  less 
than  4  x  1020  cm-3. 


To  compare  the  limitations  imposed  by  heavy  doping  on  N+P  and 
P+N  cells,  we  present,  in  Chapter  i/II,  the  results  of  a  quantitative 
parametric  study  of  the  dependence  of  VQC  in  N+P  and  P+N  silicon  solar 
cells.  These  design  parameters  considered  are  the  surface  impurity 
concentration,  the  width  of  the  quasi-neutral  emitter  region,  and 
the  emitter  surface  recombination  velocity.  We  show,  based  on  our 
computer-aided  analysis,  that  it  is  more  beneficial  to  have  a  P+  emitter 
than  an  N+  emitter  to  achieve  large  VQC  in  silicon  solar  cells. 

Finally,  we  summarize  our  studies  and  review  the  main  conclusions 
and  accomplishments  of  this  dissertation  in  Chapter  VIII.  We  also  dis- 
cuss the  scope  and  limitations  of  this  work  and  provide  suggestions 
for  future  research. 


CHAPTER  II 

AUGER  RECOMBINATION  IN  HEAVILY  DOPED  SHALLOW-EMITTER 
SILICON  PN-JUNCTION  SOLAR  CELLS,  DIODES,  AND  TRANSISTORS 


2.1  Introduction 

Because  the  recombination  current  occurring  in  the  quasi- 
neutral  emitter  can  limit  the  current  gain  of  silicon  junction  tran- 
sistors and  the  open-circuit  voltage  of  pn-junction  silicon  solar 
cells,  the  physical  origin  of  this  current  is  of  interest.  In  an  at- 
tempt to  develop  a  simple  physical  model  that  is  consistent  with 
experimental  data,  a  recent  study  [30]  asserted  that  the  inclusion  of 
Auger  recombination  in  the  heavily  doped  regions  of  the  emitter  is 
alone  sufficient  to  explain  the  data,  and  that  it  is  unnecessary  to 
include  the  heavy-doping  effect  of  bandgap  narrowing  [19].  This  con- 
clusion is  questionable  since  Auger  recombination  was  accounted  for 
qualitatively  in  the  study  [30],  rather  than  by  a  rigorous  quantitative 
treatment. 

To  remove  this  uncertainty,  we  present  in  this  chapter  a 
rigorous  analytic  evaluation  of  a  model  for  the  emitter  region  that 
includes  Auger  recombination  but  excludes  bandgap  narrowing  and 
degeneracy  of  the  charge  carriers.  We  show  that  for  silicon  pn- 
junction  solar  cells  this  model  cannot  explain  the  experimentally 
observed  values  of  VqC  [3,31]  and  the  observed  dependence  of  VQC 
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on  surface  treatment  [29].  It  follows  that  for  bipolar  transistors, 
the  model  cannot  explain  the  measured  values  of  the  common-emitter 
current  gain.  Inclusion  of  Fermi-Dirac  statistics  would  worsen  the 
discrepancy  between  theory  and  experiment.  Our  conclusions  are  con- 
trary to  the  assertions  of  [30]  but  are  consistent  with  previous 
numerical  treatments  of  the  problem  [32-35].  Our  analysis,  therefore, 
reinforces  the  importance  of  bandgap  narrowing  in  silicon  devices. 

2.2  Derivation 
We  consider  an  N-type,  heavily  doped,  inhomogeneous  emitter 
region  having  a  Gaussian  net  impurity  profile: 


NDD(x)  =  Ns  expl 


fWt 


L 


2L 


(2.1 


where  Ns  is  the  surface  impurity  concentration  and  Wc  is  the  width  of 
the  quasi-neutral  emitter  region  in  thermal  equilibrium.  The  Auger  re- 
combination lifetime  (assuming  full  ionization  of  impurity  atoms)  is 


x(x)  =  5 — — 

A     C  Nz  (x) 


(2.2) 


where  CA  is  the  Auger  coefficient  [36].  The  relationships  (2.1)  and 
(2.2)  are  used  in  the  hole  continuity  equation  (for  the  dark  condition), 


1  d  JP  (x)  =  ap(x) 
q   dx 


taI7T 


(2.3; 


where  AP  is  the  excess  hole  density.  If  we  let 


12 


y(x)  =  N  (x)  •  AP(x) 


(2.4) 


and  use  the  general  expression  for  the  hole  current  density,  including 
both  drift  and  diffusion  components,  then,  for  the  usual  low-injection 
conditions,  (2.3)  becomes 


dx2 


W£  -  x 


2LC 


WF  -  x 

dy .  exPti2l-Sr-. 

dx        I? 


y  =  o 


(2.5) 


where 


L2  a  Dp 
LA  S 


(2.e; 


is  the  square  of  the  hole  diffusion  length  at  the  emitter  surface.  In 

(2.6),  Dp  is  an  average  value  for  the  hole  diffusion  coefficient  in  the 

emitter.  The  solution  of  the  differential  equation  (2.5)  is  similar 
to  that  given  in  [37]: 


y .  A  s1nnki  erf  pfeii 


2L  j 


+  B  cosh 


!/tt  L 


err 


WE  -  X| 
2L  J 


(2.7) 


The  coefficients  A  and  B  are  determined  by  the  boundary  conditions  [11] 
2   r   .  . 


AP(0) 


Woy 


exp 


-  1 


(2.8) 


and 


Jp(WE)  =  q  •  Sp  •  AP(WE) 


(2.9) 
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In   (2.9),  Sp  is  the  hole  recombination  velocity  at  the  emitter  sur- 
face. 

The  minority-carrier  current  can  be  written  (for  low-level 
injection)  as 


i   (x)   =  -      qDP.        dy_ 


(2.10) 


Substitution  of  the  expressions  obtained  for  A  and  B  into  (2.10) 
yields  the  following  expression  for  the  emitter  saturation  current 
density: 


qDpn1 


2  cosh 


JP0 


/tt  L 


LLA 


erf 


+  (Dp/SpLA)   sinh 


n 


NCL 


SLA 


sinh 


&i  «*&!!♦ 


LA     erf|2Til  +  (VSpLA>C0Sh 


/if  L 


erf 


Wr 


V.2LJ 


{Hi  erff  E 


Equation  (2.11)   reduces  to: 


(2.11) 


P0         NSLA 


q  Dp  n,-  Vrz  , 

'     '       coth  ^-^  erf 


-A 


rwE 

2L 


(2.12) 


as  Sp  — ►  °°,  and  to  the  familiar  expression: 

2        rwr) 


q  Dp  n. 
]P0  =     NeL. 


coth 


"A 


!2.13) 


as  Sp  — *■   °°  and  L  — *■   ™   (L  — ►  °°  corresponds  to  a  flat  impurity  profile 


with  N 


DDl 


*s> 
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With  the  charge-control  relation. 


JD(0) 


Qp 


(2.14) 


we  can  evaluate  xp,  the  effective  hole  lifetime  in  the  emitter.  The 
charge  storage  in  the  emitter  is 


,WE 


Qp  =  q 


AP(x)  dx 


(2.15) 


0 


The  contribution  to  Jp(0)  due  to  surface  recombination,  JpS,  and  due 
to  Auger  recombination,  JpA,  can  be  determined  by  integrating  (2.3)  over 
the  quasi-neutral  emitter  region: 


fWE 


Jp(0)  =  Jp(W£)  +  q 


AP(x 

TA 


dx 


(2.16) 


JPS  +  JPA 


Using  the  mean-value  theorem,  we  can  write 


fwE 


^4dx  =  ^ 

TA^)     ta 


(2.17) 


Thus,  from  (2.14),  (2.16),  and  (2.17),  we  obtain 


Tp    TA    Tt 


(2.is: 
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where 


QP 
t  =tt  (2.19) 


■t  ~  J 


defines  the  hole  transit  time,  the  average  time  required  for  a  hole  to 
cross  the  quasi -neutral  emitter  region  and  recombine  at  the  surface. 
In  the  limiting  case  of  negligible  bulk  recombination  (x^  — ►  «),  tp 
is  equal  to  tt- 

To  summarize  the  analysis,  we  have  obtained  a  closed-form 
solution  of  the  continuity  equation,  and  have  derived  an  expression 
for  the  hole  current  density  when  Auger  recombination  is  the  dominant 
bulk  recombination  mechanism  in  a  heavily  doped  N-type  emitter.  Bandgap 
narrowing  and  Fermi -Dirac  statistics  have  not  been  included.  The 
surface  and  Auger  components  of  the  emitter  recombination  current  and 
the  associated  hole  transport  parameters  in  the  emitter  can  be  easily 
obtained  form  (2.14),  (2.15),  (2.16),  (2.17),  and  (2.19). 

2.3  Discussion 
The  preceding  development  is  then  an  analytic  treat- 
ment of  the  minority-carrier  transport  in  the  emitter  of  pn-junction 
devices  for  the  case  in  which  Auger  recombination  is  included  but 
bandgap  narrowing  is  neglected.  In  this  treatment  a  Gaussian  impurity 
profile  is  assumed,  and  a  finite  surface  recombination  velocity  Sp 
is  included.  Inclusion  of  a  finite  Sp  is  particularly  important  for 
solar  cells  [29], 

To  illustrate  quantitative  results  of  this  model,  we  apply 
it  to  an  N  P  silicon  solar  cell  with  a  base  doping  density 


15 


NAA  =  5  x  -I017cnf3.  We  let  N$  =  1020cnf3  and  WE  =  0.25  ym.  The 
electron  diffusion  length  l_n  in  the  base  is  assumed  to  be  *  80  um, 
consistent  with  experimental  data  [1],  In  Figure  2-1,  we  plot  JpQ 
and  its  two  components,  Jpso  and  JpA0,  as  functions  of  Sp.  For  low 
values  of  Sp  (about  103  cm/sec),  Jpso  «  JpAQ.  For  the  high  values 
of  Sp  (above  10  cm/sec)  JpsQ  *  2JpAQ.  As  Figure  2-1  shows,  the  two 
components  of  JpQ  are  equal  when  Sp  a  105  cm/sec.  When  the  emitter 
current  is  due  mainly  to  Auger  recombination  (Sp  -  103  cm/sec  or  lower), 
for  an  assumed  short-circuit-current  density  of  23  mA/cm2  [1]  ,  we 
get  V0C  ~  680  mV-  This  volta9e  disagrees  markedly  with  the  cor- 
responding maximum  observed  VgC  of  610  mV.  The  disagreement 
corresponds  to  a  discrepancy  in  JP0  of  about  a  factor  of  15. 

We  consider  now  the  variations  of  JPfj  as  a  function  of  WE, 
with  Ms  fixed  at  1020cm"3.  In  Figure  2-2,  we  plot  JpQ  versus  W£ 

2  a 

for  Sp  -  10  cm/sec  and  Sp  =  10  cm/sec.  The  figure  shows  the  dominance 
of  bulk  Auger  recombination  for  large  WE,  independent  of  the  value 
of  Sp,  and  the  dominance  of  surface  recombination  for  small  Wr  and 
large  Sp.  For  low  values  of  Sp(Sp  =  102  cm/sec  in  Figure  2-2)  the 
emitter  current  is  due  mainly  to  Auger  recombination  (Jpo  a  jpAQ)  and 
it  decreases  slightly  for  small  WE.  The  largest  value  of  the  Auger 
recombination  current  is  about  3  x  10~14A/cm2,  corresponding  to 
VgC  a  680  mV.  Again,  this  voltage  disagrees  appreciably  with  the 
610  mV  value  that  is  observed  experimentally. 

The  dominance  of  the  surface  recombination  is  emphasized 
in  Figure  2-3,  where  we  let  Sp  =  108_cm/sec  and  plot  the  resulting 


17 


o 

.(— 

u 

£5 

a; 

s 

o 

CD 

i 

0) 

o 
u 

CO 

o 

£_ 

\       «4-    <U 


o 


CL 

,r" 

0) 

CO 

•a 

c 

<T3 

03 

Q.  </> 

•"D 

a 

+J 

:> 

<1) 

o 

S- 

< 

(7ujo/v)odi? 


o 

o 
o 

o 

S- 
1 
0) 

u 

CO 

ro 
<+- 

S- 

6 

3 
a; 

N 

+j 

O 

■i—    -M 

OOJ 
CO 

tO 

n  a. 

O 

o 

-"—•• 

o 

E 

i/l    o 

LO 

=i 

S  -M 

o 

o  to 

c 

00     -i- 
O  JO 

UJ 

O 
II     c_> 

sr 

£ 

o 

O.S- 

oo    i 

E    CD 

ro 

(0    3 

< 

o 

o 

in  jz 

— »  -M 

CM 

C\J 

6 

o    o 

o 
in  d 

o 


(_ujo/v)odr 


19 


10 


-8 


10 


rS 


Sp  =  10     cm/sec 


o 
a; 

CO 


10 


rIO 


10 


J L 


J L 


0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


WE(A-'m) 


Figure  2.3     The  effectivejiole  lifetime  tr  and  its  two  components,  the 
transit  time  xt  and  the  average  Auger  lifetime  t,,  versus  Wc 
for  SP  =  1CP  cm/sec  rt  h 
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effective  hole  lifetime  and  its  components  as  functions  of  Wr.  The 

?n  -3 
surface  concentration  Ns  is  lCrucm  ..  Note  that  the  transit  time 

essentially  determines  the  emitter  recombination  current  for  WE  <0.3  um. 
This  corresponds  to  a  transparent  emitter  [24]  in  which  the  minority 
carriers  recombine  mainly  at  the  emitter  surface.  The  transparent- 
emitter  model  is  consistent  with  experimental  observations  in  con- 
ventional pn-junction  silicon  solar  cells  [29]  ,  and,  because  of 
high  recombination  velocity  at  an  ohmic  contact,  the  model  can  be 
expected  to  apply  also  to  shallow-emitter  silicon  transistors.  For 
larger  Wr,  the  average  Auger  lifetime  is  smaller  than  the  transit  time 
and  Auger  recombination  becomes  important. 

Finally,  Figure  2-4  shows  the  dependence  of  V„r  on  the  surface 

concentration  N<-  for  different  values  of  Sp  with  WE  =  0.25  ym.  When  N$ 

19   -3 
is  relatively  small  (about  10   cm  ),  surface  recombination  is  the 

dominant  recombination  mechanism  in  the  emitter,  and  Yqq  is  very 
sensitive  to  Sp.  As  N$  increases  toward  102'  cm-3,  VQC  is  eventually 
limited  by  Auger  recombination.  The  value  of  Vr~  that  results  when 
the  emitter  current  is  dominated  by  Auger  recombination  (Nc  larger 
than  2  x  10   cm)  is  about  680  mV.  Once  again,  this  result  cor- 
responds to  a  discrepancy  with  experiment  of  about  a  factor  of  15  in 
emitter  recombination  current. 

2.4  Conclusion 
The  analytic  results  we  have  presented  demonstrate  that  a 
model  for  the  emitter  of  silicon  pn-junction  devices  based  on  Auger 
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700 


Ns  (cm-3) 

Figure  2.4  Vqq  versus  the  surface  concentration  N$  for  different 
values  of  Sp  with  Wr  =  0.25  ym 
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recombination  without  bandgap  narrowing  cannot  explain  experimental 
values  of  Vg^  in  pn-junction  solar  cells.  It  follows  that  this  model 
also  cannot  explain  the  low  values  of  the  common-emitter  current  gain 
measured  in  bipolar  transistors.  In  this  demonstration,  we  have  used 
Maxwell -Boltzmann  statistics.  Use  of  Fermi-Dirac  statistics  would 
worsen  the  discrepancies  between  the  model  and  experimental  measure- 
ments [1].  Contrary  to  the  conclusions  of  [30],  therefore,  we  must 
conclude  that  physical  mechanisms,  e.g.,  bandgap  narrowing,  in  addition 
to  Auger  recombination  are  responsible  for  the  large  emitter  recombi- 
nation currents  observed  in  pn-junction  devices. 


CHAPTER  III 

RECOMBINATION  THROUGH  DEFECTS  IN  HEAVILY  DOPED 
SHALLOW-EMITTER  SILICON  PN-JUNCTION  DEVICES 


3.1  Introduction 

Having  established  in  Chapter  II  that  Auger  recombination  alone 
cannot  explain  the  observed  experimental  values  of  Vnr  in  silicon  solar 
cells,  we  must  consider  other  heavy-doping  mechanisms  that  will  bring 
theory  and  experiment  into  agreement. 

In  this  chapter,  we  consider  recombination  mechanisms,  other 
than  the  Auger  process,  that  can  possibly  occur  in  heavily  doped 
shallow-emitter  regions.  These  mechanisms  are  recombination  via  defects 
and  surface  recombination.  We  compare  the  effects  of  these  two 
mechanisms  with  ^ecent  experimental  observations  showing  the  sensitivity 
of  the  open-circuit  voltage  in  silicon  pn-junction  solar  cells  to 
surface  treatment.  These  experiments  are  discussed  briefly  in  Section 
3.2. 

In  Section  3.3,  we  discuss  Shockley-Read-Hall  recombination  via 
defects.  We  argue  that  recombination  via  defects  alone  yields  results 
that  are  inconsistent  with  experimental  observations  of  Vqq  in  silicon 
solar  cells.  If  the  Shockley-Read-Hall  lifetime  is  very  short  (in 
the  order  of  picoseconds),  it  will  yield  a  low  value  of  Vqc»  as  ob- 
served experimentally.  However,  in  that  case,  Vgr  will  be  insensitive 
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to  surface  treatment  (i.e.,  surface  recombination  velocity)  contrary 
to  the  experimental  observations  discussed  in  Section  3.2.  For  larger 
Shockley-Read-Hall  lifetimes,  which  will  yield  a  VQC  that  is  sensitive 
to  variations  in  the  surface  recombination  velocity,  the  calculated  Vqc 
is  larger  than  the  experimentally  observed  values. 

From  Section  3.3  and  Chapter  II,  we  conclude  in  Section  3.4 
that  a  short  minority-carrier  lifetime  cannot  yield  an  excessive 
recombination  current  that  agrees  with  experiment,  and  at  the  same  time 
be  consistent  with  the  recent  experimental  observations  discussed  in 
Section  3.2.  Therefore,  other  heavy-doping  mechanisms  that  can  provide 
a  large  minority-carrier  charge  storage  must  be  considered.  This 
conclusion  is  corroborated  by  recent  experimental  measurements  of  the 
short-wavelength  quantum  efficiency  of  silicon  solar  cells  [38]. 
Excessive  charge  storage  in  heavily  doped  semiconductor  regions  is 
attributed  to  either  an  effective  excess  intrinsic  carrier  density  or 
an  effective  deficit  majority-carrier  density. 

3.2  Importance  of  Surface  Recombination 
Fossum  et  al.  [29]  recently  reported  a  considerable  improvement 
in  the  open-circuit  voltage  in  a  variety  of  silicon  solar  cells  by 
reducing  the  recombination  velocity  at  the  emitter  surface.  These 
solar  calls  include  diffused,  implanted,  back-surface-field  (BSF), 
and  high-low-emitter  (HLE)  structures,  of  both  N-on-P  and  P-on-N 
types  [29]. 

It  was  observed  that  the  growth  of  a  thin  thermal  Si02  layer 
on  the  front  (emitter)  surface  increased  Vgr,  by  about  20  mV  in  some 
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cells  [29].  Table  II,  taken  from  [29],  shows  the  improvement  attained 
in  Vq-  in  the  solar  cells  described  above  by  growing  thermal  Si02 
on  the  emitter  surface. 

Other  experiments,  reported  in  [29],  indicate  that  the  thermal 
growth  of  Si02  reduced  the  surface  recombination  velocity  rather  than 
improved  the  bulk  minority-carrier  lifetime  in  the  emitter.  The  reduc- 
tion of  the  surface  recombination  velocity  by  thermal  Si02  was  observed 
for  both  iN+-  and  P+-emitters.  These  experiments  indicating  the  sensi- 
tivity of  Vqc  to  the  emitter  surface  condition  imply  that  the  emitter 
region  is  at  least  partially  transparent  to  minority  carriers.  That  is. 
a  significant  portion  of  the  emitter  recombination  current  occurs  at 
the  surface  (prior  to  the  surface  passivation). 

3.3  Shockley-Read-Hall  Recombination  via  Defects 
Defects  in  doped  silicon  can  result  from  the  diffusion  of 
impurities  and  the  subsequent  heat  treatment  of  pn- junction  devices. 
These  defects  are  due  primarily  to  vacancies  and  vacancy  complexes 
near  the  surface  region  of  the  device  [39]. 

Lindmayer  and  Allison  [40]  suggested  a  dead-layer  model  for 
the  emitter  of  thin  diffused  silicon  solar  cells  in  which  the  thin 
surface  layer  is  strongly  damaged  by  the  diffusion  of  impurities.  In 
this  surface  layer,  there  are  large  defect  densities  that  result  in 
very  short  minority-carrier  lifetimes  (in  the  picosecond  range).  The 
short  lifetimes  yield  a  large  (dark)  emitter  recombination  current 
that  limits  the  open-circuit  voltage  of  solar  cells  to  values  lower 
than  those  predicted  by  classical  pn-junction  theory. 
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In  a  treatment  of  recombination  via  defects  and  other  heavy- 
doping  effects,  Lindhoim  et  al .  [9]  studied  the  effect  on  the  dark 
emitter  current  of  process-induced  defects  in  the  emitter  of  silicon 
junction  solar  cells.  They  derived  an  impurity-concentration-dependent 
defect  density  Njj  [9]: 

NTtU)  =  K[NDD(x)  +  NM]m  (3.1) 

for  an  N  -region.  The  integer  m  can  vary  from  1  to  4  depending  on  the 
type  of  vacancy  in  the  semiconductor.  In  (3.1),  Nqq(x)  is  the 
impurity  concentration  in  an  N+-emitter,  NAA  is  the  P-type  base  impurity 
concentration,  and  K  is  a  constant.  For  an  assumed  linear  impurity 
profile  near  the  junction,  that  is,  an  impurity  concentration  that 
depends  linearly  on  distance,  they  found,  fcr  an  emitter  junction  depth 
of  0.25  ym,  that  the  excess  minority  carriers  are  packed  in  a  very  thin 

o 

layer  (about  200  A  in  width)  next  to  the  space-charge  region.  This 
situation,  illustrated  in  Figure  3.1,  means  that  most  of  the  minority 
carriers  could  not  reach  the  surface. 

The  above  conclusion  is  similar  to  earlier  bipolar- transistor 
studies  [41,42]  of  the  effect  of  a  built-in  electric  field  in  a  dif- 
fused quasi-neutral  emitter  region,  having  a  constant  minority- 
carrier  lifetime.  In  [9],  the  position-dependent  minority-carrier 
lifetime,  due  to  the  position  dependence  of  the  defect  density, 
enhanced  the  confinement  of  the  excess  injected  minority  carriers 
to  a  region  near  the  junction. 

We  consider  now  the  position-dependent  defect  density  model 
of  Lindhoim  et  al .  [9]  in  which  the  defect  density  is  proportional 
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Figure  3.1  Sketch  of  minority-carrier  concentration  in  a  heavily  doped 
emitter  having  a  graded  impurity  profile  and  a  large 
defect  density,  after  [9] 
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to  an  integer  power  of  the  doping  density  as  in  (3.1).  We  assess 
the  importance  of  defect  recombination  by  determining  if  recombination 
via  defects  is  alone  sufficient  to  explain  the  experimentally  low 
values  of  Vnp  and  its  observed  sensitivity  to  surface  treatment.  We 
consider  three  cases  for  the  integer  m  in  (3.1 ) :  m  =  1 ,  m  =  2,  and 
m  >  2. 

If  m  =  1,  the  Shockley-Read-Hall  lifetime  for  low-level 
injection  is 

1 
TSRH  "  CSRH  N(x)  (3-2) 

where  C$rw   is  a  constant.  Fossum  [43]  considered  this  case  in  a  de- 
tailed numerical  study  of  silicon  solar  cells.  He  showed  that  in  order 
to  predict  a  value  of  Vqq  comparable  to  experimentally  observed  values 
in  silicon  solar  cells,  very  short  minority-carrier  lifetimes  in  the 
emitter  region  (~  1  nsec)  are  needed.  Such  short  lifetimes  yielded 
a  Vqq  that  was  insensitive  to  variations  in  the  surface  recombination 
velocity  at  the  emitter  surface  [43],  contrary  to  experimental  obser- 
vations [29].  Therefore,  although  short  Shockley-Read-Hall  lifetimes 
that  vary  inversely  with  N(x)  can  yield  Vqc  in  the  range  of  experimental 
values,  they  cannot  account  for  the  sensitivity  of  Vqq  to  the  surface 
treatment. 

For  m  =  2,  the  Shockley-Read-Hall  lifetime  has  the  same 
dependence  on  the  doping  density  N(x)  as  for  Auger  recombination, 
except  that  the  coefficient  CSRH  is  different  than  CA: 
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1 
TSRH  "  CSRH  [N(x)]T  (3'3> 

In  this  case,  for  Shockley-Read-Hall  recombination  to  dominate  over 
Auger  recombination,  C^  must  be  larger  than  the  Auger  coefficient  Cn. 
This  yields  tsrh  <  ta> 

We  used  the  analytic  model  of  Chapter  II  to  calculate  Vqq, 
the  transit  time,  and  the  average  bulk  lifetime  for  this  case  (with 
m  =  2).  We  found  that  for  Ca  <  C<-Rm  -  4C^,  the  transit  time  is  smaller 
than  the  average  bulk  lifetime  for  high  values  of  surface  recombination 

o 

velocity  S  (about  10  cm/sec),  and  consequently  Vqq  is  sensitive  to 
variations  in  S.  However,  the  magnitude  of  Vqq  we  obtained  (>  670  mV) 
was  considerably  higher  than  the  experimentally  observed  value  of 
about  600  mV.  For  C$^  >   4CA,  our  calculations  yielded  a  bulk  lifetime 
in  the  order  of  a  fraction  of  a  nanosecond,  which  is  much  smaller  than 
the  transit  time.  Hence,  although  the  short  lifetimes  yield  lower 
values  of  Vqq  (<  670  mV),  they  result  in  a  Vqq  that  is  insensitive  to 
variations  in  S. 

For  basically  the  same  reason  discussed  above,  we  infer  that 
for  m  >  2  in  (3.1),  the  emitter  is  expected  to  be  insensitive  to  surface 
recombination.  As  the  lifetime  varies  faster  than  [N(x}]~  ,  the  injected 
minority-carrier  concentration  falls  sharply  with  increasing  distance 
away  from  the  junction.  Thus,  most  of  the  injected  minority  carriers 
cannot  reach  the  emitter  surface  before  they  recombine  in  the  bulk 
via  defect  centers. 
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3.4  The  Excess  Minority-Carrier 
Charge  Storage  in  Heavily  Doped  Regions 

We  discussed  Shockley-Read-Hall  recombination  via  defects 
in  the  last  section  and  showed  that  for  Shockley-Read-Hall  lifetimes 
that  have  an  inverse-power  functional  dependence  on  the  doping  density, 
recombination  via  defects  cannot  alone  be  the  dominant  mechanism 
limiting  Vqq  in  silicon  solar  cells.  This  is  because  such  a  recom- 
bination mechanism  cannot  yield  results  that  are  simultaneously  in 
agreement  with  the  low  values  of  Vqq  observed  in  silicon  solar  cells 
and  with  the  sensitivity  of  Vgr  to  surface  treatment.  Furthermore, 
recent  experimental  investigation  of  the  short  wavelength  spectral 
response  of  silicon  solar  cells  indicates  that  heavy-doping  effects, 
other  than  short  lifetimes,  are  required  to  obtain  agreement  between 
experimental  and  theoretical  results  [38]. 

Combining  the  above  conclusion  with  the  conclusion  of  Chapter 
II,  we  deduce  that  short  minority-carrier  lifetimes  cannot  be  the  only 
cause  of  the  large  emitter  current  in  silicon  solar  cells.  A  similar 
conclusion  is  anticipated  for  bipolar  transistors.  From  Equation  (1.1), 
we  infer  that  the  minority-carrier  charge  storage  in  heavily  doped 
regions  must  then  be  large  to  account  for  the  large  current. 

We  now  show  qualitatively  that  a  large  minority-carrier  charge 
storage  is  due  either  to  an  excess  intrinsic  carrier  density  or  to  a 
deficit  majority-carrier  concentration. 

The  minority-carrier  charge  storage  Qp,  in  an  N-type  region 
for  example,  is 
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Qn  =  q 


lP(x)  dx 


:3.4) 


0 

where  AP(x)  is  the  excess  minority-carrier  (hole)  concentration.  For 
non-equilibrium  conditions  in  a  non-degenerate  semiconductor, 


PN  =  n.j  exp 


EFN  "  EFP 


L 


kT 


where  Epp  and  Epj|  are  the  hole  and  electron  quasi -Fermi  levels, 
From  (3.5),  the  condition  of  low-injection, 


(3.5) 


N  *  N 


0 


(3.6) 


and  the  fact  that 


P  -   AP 


(3.7) 


we  get 


AP  -  M^ 

N0 


(3.8) 


If  we  assume  full  ionization  of  impurities, 


%   *  NDD 


(3.9) 


and  (3.8)  becomes 

2 


AP 


DD 


(3.10: 


32 


If,  because  of  some  physical  mechanism,  e.g.,  energy-bandgap  narrowing, 
n^r  is  larger  than  its  intrinsic  value  (in  pure  silicon),  or  N«  is 
smaller  than  the  doping  concentration  because  of,  e.g.,  incomplete 
ionization  of  the  impurity  atoms,  then  AP,  and  consequently  Qp  in  (3.4), 
will  be  larger  than  anticipated,  and  thus,  by  (1.1),  the  emitter  current 

will  be  excessively  large. 

2 

The  two  effects  mentioned  above,  n^  larger  than  its  intrinsic 

value  and  Nq  <  N™,  can  be  phenomenologically  incorporated  into  an 
effective  excess  intrinsic  carrier  density  n?  such  that  (3.10)  can  be 
written  as 


2 

AP  ccjjli  (3.11) 

NDD 

for  a  non-degenerate  N-type  semiconductor.  For  a  degenerate  semiconductor, 

we  show  later  that  n?  can  include  the  effects  of  Fermi-Dirac  statistics, 
ie 

and  hence  (3.11)  is  still  applicable.  Note  that,  in  this  way,  heavy- 

o 

doping  effects  can  be  incorporated  in  a  single  parameter,  n*-  . 

In  the  next  chapter,  we  will  consider  deionization  of  impurities, 

as  discussed  by  the  recent  papers  of  Heasell  [25,44]  and  Popovic  [45], 

2    2 
as  a  possible  mechanism  that  results  in  n^  >  n^. 
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CHAPTER  IV 

ON  THE  DEIONIZATION  OF  IMPURITIES  AS  AN  EXPLANATION  FOR 
EXCESS  INTRINSIC  CARRIER  DENSITY  IN  HEAVILY  DOPED  SILICON 


Experimental  measurements  of  the  collector  currents  in  silicon 
bipolar  transistors  [20]  have  indicated  that  n-  ,  the  square  of  the 

effective  intrinsic  carrier  concentration,  increases  to  values  well 

2 
above  n.,  the  classical  intrinsic  carrier  concentration,  as  the  doping 

level  in  the  quasi-neutral  base  region  increases.  This  was  interpreted 

by  Slotboom  and  DeGraaf f  to  be  due  to  energy-bandgap  narrowing  in  the 

quasi-neutral  base  region,  which  was  heavily  doped  in  the  devices  used 

for  their  study  [20]. 

Recently,  Popovic  [45]  presented  an  alternate  interpretation 

2 
for  the  increase  of  n^  with  the  doping  level.  He  argued  that,  in 

the  range  of  impurity  concentrations  of  10'  cm   to  10   cm   in 

2 
silicon,  the  increase  of  n^  can  be  attributed  to  a  Fermi  energy 

difference  AE  corresponding  to  different  density-of-states  models, 

whose  physical  significance  will  be  discussed  later.  In  recent  papers  [25, 

44],  Heasell  gave  another  interpretation  based  on  the  assumption  of 

a  constant  impurity  activation  energy  at  high  doping  levels.  The 

models  of  Popovic  and  Heasell  resemble  each  other  in  that  they  both 

emphasize  the  importance  of  deionization  of  impurity  levels,  which 

results  in  the  majority  carrier  density  being  substantially  lower  than 

the  net  doping  density. 
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We  discuss  here  in  detail  the  models  of  Popovic  and  Heasell, 
and  show  that  their  underlying  assumptions  are  inconsistent  with  a 
variety  of  experimental  observations  in  heavily  doped  silicon.  Hence 
we  conclude  that  these  interpretations  do  not  dislodge  energy-bandgap 
narrowing  as  a  possible  mechanism  underlying  the  increase  in  n^e  in 
heavily  doped  silicon. 

As  the  basis  for  his  model,  Heasell  assumed  (a)  that  the 
bandgap  of  silicon  is  not  affected  by  heavy  doping,  and  (b)  that  the 
impurity  energy-levels  in  the  bandgap  due  to  the  presence  of  shallow- 
level  impurity  atoms  &re   independent  of  the  impurity  concentration. 
For  heavily  doped  silicon,  these  assumptions  lead  to  a  small  ionization 
ratio,  which  we  define  as  the  ratio  of  the  concentration  of  ionized 
impurity  atoms  to  the  total  concentration  of  impurity  atoms.  To  show 
this,  consider  a  donor  energy  level  of  En.  Then  the  majority-carrier 
concentration  N  (electrons  for  an  N  region),  which  equals  the  ionized 
impurity  concentration,  is  related  to  the  total  donor  concentration, 
Nqq,  for  uncompensated  N+  silicon  by  £46] 

N  -  [NDDNc/gD]1/2exp[-(Ec  -  ED)/2kT]  (4.1) 

18    3 

for  NqD  >  10   cm  .  Here  gn  is  the  impurity  energy-level  degeneracy 

factor,  Np  is  the  conduction-band  effective  density  of  states,  and 
Ef  is  the  conduction-band  edge.  Thus,  in  the  Heasell  model,  where  Eq 
does  not  depend  on  the  impurity  concentration,  the  majority-carrier 
concentration  varies  in  proportion  as  the  square  root  of  the  impurity 
concentration.  If  Nnn  increases  by  two  orders  of  magnitude,  N  increases 
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by  only  one  order  of  magnitude.  Hence,  the  NP  (electron-hole)  product 
(given,  at  equilibrium  for  low  doping  concentrations,  by  the  law  of 
mass  action)  implies  a  larger  minority-carrier  density  than  is 

predicted  classically.  This  has  the  same  effect  in  device  analysis 

2  2 

as  would  n.j  being  larger  than  n^. 

The  above  assumption  of  a  fixed  discrete  impurity  energy  level 

disagrees  with  the  experimental  observations  of  Pearson  and  Bardeen 

[47]  and  Penin  et  al.  [48].  They  found  that  the  ionization  energy 

of  the  impurity  atoms  vanishes  as  the  impurity  concentration  approaches 

18-3  - 

about  3  x  10   cm   fcr  both  n-type  and  p-type  silicon.  Therefore, 

Heasell's  assumption  of  a  fixed  impurity  energy  level  in  the  bandgap 

of  heavily  doped  silicon  conflicts  with  long  established  experimental 

findings. 

To  modify  Heasell's  approach  to  bring  it  into  harmony  with 

the  results  of  these  experiments,  we  have  calculated  the  ionization 

ratio  for  heavily  doped  silicon  using  an  impurity-concentration 

dependence  of  the  ionization  energy  based  on  the  experimental  results 

of  Pearson  and  Bardeen  [47].  We  found,  in  agreement  with  [49],  that 

the  lowest  ionization  ratio  at  room  temperature  is  about  90%  for 

phosphorus-doped  silicon,  and  about  70%  for  boron-doped  silicon. 

This  difference  in  the  degree  of  ionization  arises  from  the  difference 

in  the  effective  density  of  states  in  the  conduction  and  valance  bands 

and  in  the  degeneracy  factors  applying  to  the  donor  and  acceptor 

impurity  energy  levels.  This  deionization  of  impurities,  based  on  the 

experimental  results  of  Pearson  and  Bardeen  [47],  is  not  large  enough 


2 
to  account  for  the  experimentally  observed  values  of  n.  quoted  by 

le 

Heasell  [25].  Thus  Heasell's  assertion  that  deionization  is  the 
first-order  mechanism  underlying  the  large  values  of  n?  observed  in 
heavily  doped  silicon  is  not  justified. 

We  discuss  now  the  treatment  of  Popovic  [45].  Popovic 
related  the  increase  of  n?  observed  in  the  base  region  of  a  bipolar 
transistor  as  the  impurity  concentration  increases  above  10^  cm 
[20]  to  a  difference  in  Fermi  energy  AE.  The  energy  difference  AE 
originates  from  calculations  based  on  two  different  models  for  the 
density  of  states  in  heavily  doped  silicon. 

The  first  (classical)  model  assumes  parabolic  conduction  and 
valance  bands;  i.e.,  the  density  of  states  is  assumed  proportional  to 
the  square  root  of  energy,  and  full  impurity  ionization  is  assumed. 
This  model  yields  a  Fermi  energy  Ep,  calculated  from  the  quasi- 
neutrality  condition,  that  increases  monotonically  as  the  doping  level 
increases.  Above  an  impurity  concentration  of  about  10   cm"  ,  Ep 
enters  the  conduction  band  (for  n-type  material)  or  the  valance  band 
(for  p-type  material). 

The  second  model  assumes  that  (a)  the  conduction-band  and 
valance-band  density  of  states  have  exponential  tails  that  extend 
into  the  forbidden  bandgap,  (b)  the  impurity  density  of  states  has  a 
Gaussian  form,  and  (c)  the  total  density  of  states  is  the  envelope 
of  the  density  of  states  in  (a)  and  (b)  [16].  When  used  with  the 
quasi-neutral ity  condition  in  the  base  region,  this  model  predicts 
a  different  Fermi  energy  Ep  than  that  (Ep)  calculated  from  the  classical 
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model.  The  Fermi  energy  Ep  saturates  for  high  impurity  concentrations 
at  the  low-impurity-concentration  value  of  the  impurity  energy  level 
in  the  bandgap  [16]. 

Popovic  claims  that  the  energy  difference, 

AE  =  EF  -  Ep  (4.2) 

is  an  alternate  interpretation  to  energy-bandgap  narrowing  in  silicon 
in  the  doping  range  of  1017  to  10   cm-3,  an  interpretation  that  can 
account  for  the  increase  in  n^  in  the  same  way  as  bandgap  narrowing, 
i.e., 

n?e  =  n^  exp(AEAT)  (4.3) 

for  Maxwell -Boltzmann  statistics.  The  energy  difference  AE  increases 
as  the  impurity  concentration  increases  because,  in  the  second  model 
[16]  for  the  density  of  states,  EA  saturates  at  the  impurity  energy 
level ,  while,  according  to  the  classical  (first)  model,  Ep  increases 

monotonically  as  the  impurity  concentration  increases. 

2 
To  estimate  quantitatively  the  effect  of  AE  on  n.  ,  we  con- 

1  q   _3 
sider  a  concentration  of  10   cm   phosphorus  atoms  in  silicon  at  room 

temperature.  Then,  AE  -  0.045  eV  since  Ep  is  yery   close  to  the  con- 
duction band  edge,  while  Ep  is  nearly  at  the  impurity  energy  level 
of  phosphorus  in  silicon--0.045  eV  below  the  conduction  band  edge. 
Using  this  value  of  AE,  we  obtain 

n?e  »  6n?  (4.4) 
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Since  0.045  eV  is  the  maximum  value  of  AE  in  the  range  of  doping  con- 

1 7    °     1 9   -3  ? 

sidered  (from  10   cm"J  to  10   cm  ),  ine  maximum  value  of  x\.     due  to 

the  difference  in  Fermi  energy  is  less  than  an  order  of  magnitude 
larger  than  n?.  If  we  use  this  value  of  n?  ,  for  the  range  of  doping 
specified  above,  add  the  effect  of  the  bandgap  narrowing  model  of 
Vol 'fson  and  Subashiev  [50]  (for  NDD  >  2  x  1019  cm"3),  as  suggested  by 
Popovic,  and  include  Fermi -Dirac  statistics,  we  get  values  of  the 
open-circuit  voltage  in  silicon  solar  cells  and  the  common-emitter 
gain  in  silicon  bipolar  transistors  that  are  much  higher  than  experi- 
mentally  observed  values.  Thus  the  increase  of  n^  due  to  the  satura- 
tion of  El,  assuming,  in  this  part  of  the  discussion,  that  it  saturates, 
is  too  small  to  account  for  these  experimental  findings  in  silicon 
devices. 

We  now  turn  to  the  basic  assumptions  of  Popovic 's  treatment: 
the  saturation  of  the  Fermi  level  Er  and  the  resulting  tacit  assump- 
tion of  impurity  deionization  in  heavily  doped  silicon. 

From  a  theoretical  viewpoint,  saturation  of  the  Fermi  level 
in  the  bandgap  is  due  primarily  to  impurity  band  widening,  that  is,  to 
the  spread  of  impurity  energy  levels  from  discrete  to  quasi-continuous 
levels  in  the  bandgap  at  high  impurity  concentrations  (in  the  range 
of  10   to  10   cm-3).  Impurity  band  widening  introduces  a  large 
number  of  energy  levels  in  the  bandgap  that  can  be  populated  by  the 
majority  carriers.  This  lessens  the  density  of  carriers  in  the 
majority  allowed  band,  and  makes  the  Fermi  level  tend  to  saturate. 

Deionization  of  impurities  is  tacitly  assumed  in  Popovic's 
treatment  because  the  charge  carriers  populating  the  impurity  band 
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are  regarded  to  be  localized  spatially.  This  is  so  because,  if  the 
carriers  occupying  the  impurity  band  were  not  localized,  i.e.,  if  their 
wavefunctions  were  to  extend  all  over  the  crystal,  there  would  be, 
in  effect,  energy-bandgap  narrowing  since  allowed  band-like  states 
would  exist  at  lower  energies  (for  N-type  material)  than  those  present 
in  pure  silicon.  However,  since  Popovic  asserts  in  his  model  that 

the  bandgap  does  not  change  because  of  the  presence  of  impurities 

17     19   -3 
in  the  doping  range  of  10   to  10   cm  ,  Popovic 's  model  must  assume 

that  the  impurity  band  states  are  localized  spatially.  Such  localiza- 
tion results  in  a  majority  carrier  concentration  N  that  is  considerably 
less  than  the  impurity  concentration  Nqq.  Thus,  popovic 's  model  leads 
to  the  deionization  of  impurity  atoms,  which  was  a  central  physical 
mechanism  also  in  the  Heasell  model.  In  the  Popovic  model,  impurity 
deionization  and  Fermi-level  saturation  go  hand-in-hand. 

But  we  now  note  that  Fermi  level  saturation  is  not  the  trend 
observed  experimentally  in  heavily  doped  silicon.  A  variety  of  ex- 
perimental observations  in  heavily  doped,  single-crystal  silicon 
indicates  that  the  localization  of  charge  carriers  is  greatly  reduced 
by  the  increase  in  impurity  concentration  and  the  electrons  (or  holes) 
show,  experimentally,  physical  properties  that  are  in  good  agreement 
with  the  degenerate  free-electron  gas  model  in  metals.  These  experi- 
mental observations  in  heavily  doped  silicon  involve  measurements  of: 
resistivity  and  Hall  effect  [47-48,  51-55],  electron  spin  resonance 
(ESR)  [53],  magneto-resistance  [54],  electronic  specific  heat  [56], 
and  nuclear  magnetic  resonance  (NMR)  [57-53].  A  thorough  discussion 
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of  these  mechanisms  in  group  IV  semi-conductors  is  given  in  [58].  We 
consider  briefly  here  the  implication  of  measurements  of  the  electronic 
specific  heat  [56]  and  nuclear  magnetic  resonance  [57-58]  as  they 
relate  to  the  position  of  the  Fermi  level  in  heavily  doped  silicon. 
Measurements  of  the  electronic  specific  heat  in  phosphorus- 
doped  silicon  (Nqq  >  6  x  10   cm)  show  that  the  electrons  behave  as 
a  degenerate  electron  gas  [56].  The  electronic  specific  heat  of  the 
heavily  doped  silicon  samples  (a)  has  a  linear  dependence  on  tempera- 
ture, (b)  increases  considerably  as  the  impurity  concentration 

1/3 
increases,  and  (c)  has  an  N™  dependence  on  the  impurity  concentration 

in  agreement  with  the  dependence  predicted  by  the  degenerate  free- 
electron  gas  model  for  metals  [59].  The  above  results  indicate  that, 
for  Nqq  >  6  x  10  ,  a  parabolic  density  of  states  in  the  conduction 
band  is  probably  a  good  approximation.  Furthermore,  although  no 
quantitative  estimates  of  the  degree  of  ionization  can  be  concluded  from 
measurements  of  the  electronic  specific  heat,  these  measurements  do  in- 
dicate that  almost  all  donor  electrons  occupy  delocalized  states; 
otherwise  the  specific  heat  would  not  exhibit  the  metallic  character 
described  above. 

Experimental  determination  of  the  position  of  the  Fermi  level 
and  of  the  derealization  of  donor  electrons  (that  is,  of  the  ioniza- 
tion of  impurity  atoms)  is  provided  by  nuclear  magnetic  resonance  (NMR) 
experiments  in  phosphorus-doped  silicon  [57-58].  These  experiments 
indicate,  via  the  spin  interactions  of  electrons  and  impurity  nuclei, 

1  Pi  3 

that,  at  an  impurity  concentration  of  about  3  x  10   cm"  ,  electrons  in 
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heavily  doped  silicon  are  delocalized,  in  agreement  with  the  results 
of  Pearson  and  Bardeen  [471.  Furthermore,  at  a  phosphorus  concentration 
of  about  2  x  10 Iy  cm"  ,  the  Fermi  level  in  heavily  doped  silicon  is  at 
the  conduction  band  edge  of  silicon  [58],  that  is,  the  Fermi  level 
enters  the  conduction  band.  Thus  NMR  experiments  indicate  that  the 

Fermi  level  does  not  saturate  in  the  bandgap  of  silicon  [58].  They 

19   -3 
indicate  also  that  at  impurity  concentrations  above  2  x  10   cm 

the  impurity  band  merges  almost  completely  with  the  conduction  band 

[58],  resulting  in  a  composite  density  of  states  having,  to  a  good 

■J  10 

approximation,  a  parabolic  dependence  on  energy  (E  '  )  [57]. 

Therefore,  we  conclude  that  the  basic  assumptions  of  Popovic's 
treatment  are  not  justified  since  experiments  predict  that  the  Fermi 
level  does  not  saturate,  and  impurity  deionization is  not  a  first-order 
mechanism  in  heavily  doped  silicon. 

Our  approach  has  been  to  dispute  the  conclusions  of  Heasell 
and  Popovic  because  their  initial  premises  conflict  with  experimental 
observations.  In  the  case  of  Heasell 's  model,  we  also  showed  that  a 

version  of  it,  modified  to  be  consistent  with  experiment,  still  failed 

2 
to  explain  the  observed  variation  of  n^e  with  doping  level.  In  this 

chapter  we  have  not  discussed  theoretical  grounds  for  the  model 

o 

advanced  earlier  to  explain  the  variation  of  n^  ,  that  is,  the  model 
of  energy-bandgap  narrowing.  Our  position  is  that  bandgap  narrowing 
is  a  possible  valid  mechanism  influencing  n.  --and  thus  the  common- 
emitter  current  gain  of  silicon  transistors  and  the  open-circuit 
voltage  of  silicon  solar  cells.  But  we  believe  that  the  stage  now 
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reached  in  theoretical  and  experimental  work  on  bandgap  narrowing  does 
not  allow  a  consensus  of  agreement  to  be  reached  that  it  is  the  dominant 
mechanism  underlying  the  observed  variations  in  n?e. 

In  accordance  with  the  above  comment,  the  expression  "bandgap 
narrowing"  will  mean,  in  the  following  chapters,  a phencmenological 
narrowing  of  the  bandgap  characterized  by  an  excess  intrinsic  carrier 
concentration  nje  >  nj.  That  is,  our  use  of  the  term  "bandgap 
narrowing"  does  not  necessarily  refer  to  the  true  shrinkage  in  the 
bandgap  of  the  semiconductor.  As  shown  in  Chapter  II  and  later  in 
this  dissertation,  use  of  a  phenomenclogical  bandgap  narrowing  in  the 
analysis  of  heavily  doped  silicon  devices  is  important  to  bring  theory 
and  experiment  into  agreement. 

In  the  next  chapter,  we  incorporate  bandgap  narrowing  into 
the  analysis  of  heavily  doped  emitters  of  silicon  pn-junction  devices. 


CHAPTER  V 

HEAVILY  DOPED  TRANSPARENT-EMITTER  REGIONS 
IN  SILICON  JUNCTION  SOLAR  CELLS,  DIODES,  AND  TRANSISTORS 


5.1  Introduction 

Excess  minority  carriers  injected  into  the  emitter  of  pn- 
junction  devices  reccmbine  in  the  bulk  and  at  the  surface  of  the 
emitter.  If  the  emitter  junction  is  shallow  enough,  the  minority 
carriers  can  cross  the  quasi-neutral  emitter  region  without  appreciable 
bulk  ecombination.  The  minority  carriers  then  recombine  at  the 
emitter  surface.  For  this  case  the  emitter  is  transparent  to  the 
injected  minority  carriers,  and  an  important  parameter  then  is  the 
surface  recombination  velocity  S  at  the  emitter  surface. 

This  parameter  is  particularly  important  for  pn-junction 
silicon  solar  cells  in  which  most  of  the  illuminated  surface  is  not 
covered  by  metal.  In  devices  in  which  thermal  S i O2  covers  this 
nonmetallized  portion  of  the  surface,  experiments  show  that  S  can  be 
less  than  104  cm/sec  for  both  P-cells  (P-type  substrate)  [29]  and 
N-cells  [29].  This  value  of  S  is  orders  of  magnitude  less  than  that 
of  an  ohmic  contact  and  is  consistent  with  values  determined  earlier 
by  different  experimental  methods  [60].  Furthermore,  recent  experi- 
ments involving  P-on-N  and  N-on-P  cells  without  thermal  Si02  demon- 
strate that  the  emitter  can  be  completely  transparent  [29]. 
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The  purpose  of  this  chapter  is  to  provide  an  analytical 
treatment  of  transparent  emitter  devices,  particularly  solar  cells, 
that  is  more  complete  than  treatments  previously  available 
In  this  treatment,  we  include  the  effects  of  (a)  bandgap  narrowing 
[19,61],  (b)  Fermi-Dirac  statistics,  (c)  built-in  field  due  to  the 
impurity  profile,  and  (d)  a  finite  surface  recombination  velocity  S. 
Detailed  numerical  studies  including  these  various  effects  have  been 
done  [32,62-63],  but  they  have  not  treated  the  case  of  the  transparent 
emitter. 

A  major  result  of  this  chapter  is  the  demonstration  that  the 
transparent-emitter  model  can  predict  experimental  values  of  Vqq 
observed  on  N+P  thin  diffused  junction  silicon  solar  cells  made  on 
low-resistivity  (0.1  ft-cm)  substrates.  Thus,  the  transparent-emitter 
model  is  shown  to  provide  an  explanation  for  the  discrepancy  between 
the  prediction  of  simple  classical  theory  (Vqq  =  700  mV)  and  the 
measured  maximum  value  (Vqq  =  600  mV).  The  transparent-emitter  model 
gives  VqC  =  600  mV  for  high  values  of  Sp  (Sp  >  10  cm/sec)  provided 
the  effects  of  bandgap  narrowing  (modified  by  Fermi-Dirac  statistics) 
are  included.  This  result  suggests  that  VqC  can  be  increased  toward 
the  classical  value  of  700  mV  if  Sp  is  decreased  and  the  effects  of 
bandgap  narrowing  are  reduced.  This  is  accomplished  in  the  HLE  solar 
cells,  early  verions  of  which  have  shown  increases  in  VqC  to  the 
640-650  mV  range  [64-65]. 

In  addition  to  the  development  of  the  theory  for  the 
transparent-emitter  device,  and  its  application  to  solar  cells,  this 
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chapter  will  include  a  test  for  the  self-consistent  validity  of  the 
transparent-emitter  model.  This  test  compares  the  calculated  transit 
time  of  minority  carriers  across  the  emitter  with  the  Auger-impact 
minority-carrier  lifetime  within  the  emitter  region. 

5.2  Derivation 
We  consider  an  N-type  heavily  doped  quasi-neutral  emitter 
region;  analogous  results  apply  to  P-type  emitters.  The  minority 
carrier  current  density  in  the  N-type  region  is 

Jp(x)  =  qypAP(x)  E(x)  -  qDp^^-  (5.1) 

in  which  E(x)  is  the  thermal  equilibrium  value  of  the  effective  field 
which,  for  low-level  injection,  is  given  by 

Dp    i     dP  (x) 
^^•yxT'-^x-  ^ 

We  now  define  an  effective  intrinsic  density  n-  such  that 

n*e(x)  =  PQ(x)  NQ(x)  (5.3) 

in  which  Pq(x)  and  NQ(x)  are  the  hole  and  electron  concentrations  in 
thermal  equilibrium.  The  parameter  n?  depends  on  position  for  two 
reasons: 

1.  The  influence  of  Fermi-Dirac  statistics,  and 

2.  The  influence  of  bandgap  narrowing 

These  influences  are  discussed  in  Section  5.3.  For  Maxwell - 

2 
Boltzmann  statistics  and  no  bandgap  narrowing,  n.  is  the  square  of 
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the  intrinsic  carrier  concentration  in  silicon  and  is  a  function  of 
temperature  only. 

If  the  expressions  is  (5.2)  and  (5.3)  are  used  in  (5.1),  we 
get,  after  some  manipulations, 

r 


Nn(x) 
Jp(x)  -y dx  =  -qDpd 

n1e(x) 


ap(x; 


N0(x)" 


L 


n.  (x 
ie  ■ 


(5.4) 


If  we  integrate  (5.4)  over  the  quasi-neutral  emitter  region,  we  get 

,Wr 


ft 


N0(x 


r 


o 


Jp(x)  •  -j dx  =  -qDp 

nie(x) 


AP(x) 


Nn(x)l 


nie(x) 


(5.5) 


-0 


where  Dp  is  some  average  value  of  Dp.  If  the  emitter  is  transparent 
(transit-time  limtied),  that  is,  if  the  minority  carrier  transit  time 
Tt  is  much  less  than  the  minority  carrier  lifetime  x  (for  an  M-type 
emitter),  then  Jp  is  constant  independent  of  position  in  the  emitter. 
Use  of  the  minority  carrier  boundary  conditions  [11] 


AP(0)  =  P0(0)(exp(qV/kT)  -1) 

at  the  edge  of  the  emitter  space  charge  region,  and  [11] 

Jp(WE)  =  q  •  Sp  •  AP  (We) 

at  the  emitter  surface,  enable  (5.5)  to  be  expressed  as 
qDp  (exp(qV/kT)-  1) 


(5.e; 


(5.7) 


Jd  = 


WE  NQ(x)       Dp  ;  NQ(WE) 
-o dx  +  x 


(5.s; 


0J   nie(x)     Sp  •  n^(WE) 
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Equation  (5.8)  is  the  general  experssion  for  the  minority-carrier 
current  in  a  transparent  emitter. 

To  check  the  condition,  x  «  x  ,  required  for  transparancy, 
we  must  determine  the  steady-state  transit  time  xt,  which  is  defined 
by  the  charge  control  relation, 


t  "  Jr 


(5.9; 


Here  Qp  is  the  excess  minority  carrier  charge  storage  in  the  emitter: 


rwE 


QP  =  q 


AP(x)  dx 


15.10) 


Using   (5.4),    (5.8),   and   (5.10)   to  express  Q   ,   and  combining  with  (5.9), 
we  find  the  following  expression  for  the  minority  carrier  transit 
time: 

'-     fwE 

N0(W£)_      ^    I 
J 


fWE 


Nq(x) 


dx  + 


SP  '  nfe<W  ) 


L     0 


n2  (x; 

ie 
l^TxT 


dx 


J      0 


Wr 

f"E    rx 


0       0 


Nn(x')        nf   (x) 
0  ie  dx'   dx   > 


Nnu; 


n2   (x1 
ie 


(5.11 
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Some  special  cases  are  of  interest.  For  a  flat  impurity 
concentration  profile,  the  above  expression  reduces  to 


4    w 

-t  =  ii+v  (5-12: 

If,  furthermore,  Sp  is  infinite,  (5.12)  reduces  to  the  familiar  ex- 
pression, 


4 

5.3  Heavy-Doping  Effects 

In  thermal  equilibrium,  heavy-doping  concentrations  of  shallow 
level  impurities  affect  the  minority-carrier  concentration  in  a  quasi  - 
neutral  region  by  two  mechanisms:  bandgap  narrowing  and  Fermi -Dirac 
statistics.  These  two  mechanisms  affect  the  minority-carrier 
concentration  in  opposite  ways.  For  any  given  position  of  the  Fermi 
level  relative  to  the  band  edges,  bandgap  narrowing  tends  to  increase 
the  minority-carrier  concentration,  while  inclusion  of  Fermi-Dirac 
statistics  tends  to  decrease  the  minority-carrier  concentration  below 
the  value  calculated  using  Maxwell -Boltzmann  statistics.  The  dominance 
of  either  of  the  two  effects,  at  any  specific  impurity  concentration, 
depends  en  the  model  of  bandgap  narrowing  adopted. 

In  this  treatment,  we  assume  that  bandgap  narrowing  occurs 
without  changing  the  parabolic  dependence  on  energy  of  the  density  of 
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states  in  the  conduction  and  valance  bands.  This  is  the  rigid-band 
approximation;  it  is  discussed  in  Appendix  F. 

The  effects  of  bandgap  narrowing  and  Fermi-Dirac  statistics 
can  be  lumped  into  a  position-dependent  effective  intrinsic  carrier 
concentration  at  thermal  equilibrium  given  by 

2  ,  .   n?  exp(AEp(x)/kT) 
nie(x)=  !  +  C(n)  exptnF  (5J4) 

where 

n  =  nc  =  (Ep  -  Ec)/kT  (5.15) 

for  N-type  material,  and 

n.  =  nv  =  (Ev  -  Ep)/kT  (5.16) 

for  P-type  material.  The  derivation  of  (5.14)  appears  in  Appendix  C. 
In  (5.15)  and  (5.16),  Ec  and  Ey  are  the  edges  of  the  conduction  and 
valence  bands,  respectively,  and  Ep  is  the  Fermi  level.  The  factor 
C(n)  is  a  function  of  t),   which,  for  n  5  4  (e.g.,  NQ  <  2  x  1020  cm-3  in 
N-type  silicon) ,  is  [24] 

C(n)  =  -0.04n  +  0.3  (5.17) 

The  above  approximation  of  C(n)  gives  values  of  the  Fermi-Dirac 
integral  of  order  1/2  with  less  than  4%  error. 

In  nonequilibrium  conditions,  bandgap  narrowing  increases  the 
minority-carrier  current  by: 
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1.  Increasing  the  minority-carrier  concentration 

2.  Decreasing  the  retarding  built-in  electric  field 
acting  on  the  minority  carriers 

The  increase  in  the  minority-carrier  concentration  P  results 

o 

from  the  increase  in  n^  (x).  The  decrease  of  the  built-in  electric 

field  is  due  to  the  position  dependence  of  n.  (x)  (and  hence  of  the 
effective  bandgap)  in  the  inhomogeneously  doped  emitter. 

To  develop  a  simple  expression  illustrating  the  reduction  of 
the  electric  field,  we  now  include  only  bandgap  narrowing  excluding  the 
effect  of  Fermi-Oirac  statistics  for  the  present.  Then  the  effective 
electric  field  acting  on  the  minority  carriers,  given  in  (5.2),  can  be 
expressed  by  using  (5.3)  and  (5.14)  as: 


Dr 


£W  =  R(N0M  -  ^  •  N^ 


dN0W 


dx 


(5.13) 


where 


R(Nn) 


1  - 


N0(x)   dAEG(xT 


kT 


dNQ(x) 


(5.19! 


The  factor  R(N  )  measures  the  reduction  of  the  built-in  electric  field 
due  to  bandgap  narrowing.  For  any  model  of  bandgap  narrowing,  R(NQ) 
is  always  less  than  one.  Figure  5.1  shows  R(N  )  as  a  function  of  the 
electron  (majority  carrier)  concentration  for  three  models  of  bandgap 
narrowing:  Slotboom  and  DeGraaff  [20],  Mauser  [12],  and  Lanyon  and 
Tuft  [13]. 
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In  the  absence  of  bandgap  narrowing,  the  holas  experience 
a  retarding  electric  field  i.i  an  N-type  diffused  emitter.  liquations 
(5.18)  and  (5.19)  indicate  that  the  position  dependence  of  the  bandgap 
narrowing,  in  effect,  decrease?  the  retarding  electric  field.  The 
more  it  is  decreased  the  smaller  is  the  transit  time  for  a  specific 
surface  recombination  velocity.  In  (5.11)  the  transit  time  is  shown 
to  be  a  function  of  n^  .  In  Figure  5.2  the  transit  time  is  plotted 
as  a  function  of  W.-,  the  width  of  the  quasi-neutral  emitter  region, 
in  two  cases:  neglecting  bandgap  narrowing,  and  including  bandgap 
narrowing  (Slotboom  and  De  Graaff  model).  Mote  that  inclusion  cf 
bandgap  narrowing  makes  the  transit  time  close  tc  the  value  it  has 
if  the  impurity  profile  is  flat.  In  general,  bandgap  narrowing 
decreases  the  transit  time  if  the  impurity  profile  is  not  flat. 

Conversely,  inclusion  of  Fermi-Dirac  statistics  increases  the 
transit  time  as  can  be  seen  in  Figure  5.2  Inclusion  of  Kermi-Dirac 
statistics  shifts  the  value  of  the  transit  time  closer  to  that  calcu- 
lated when  heavy-doping  effects  are  neglected. 

5.4  Discussion 
From  (5.8),  the  minority-carrier  saturation  current  for  a 
transparent  emitter  is 

JP0  " 15.20 

w 

E  N0(x)   .  .  DpNn(Wp) 

OX  T  _  \        ' 


J   nie<¥     Spnie(WE) 


0.2 


0.4 
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0.8 


VL  (/*m) 


£ 

Figure  5.2  The  transit  time  rt  versus  the  width  of  the  emitter  region 
Wr  for  S  =  5  x  105  cm/sec  and  a  Gaussian  profile  with:  no 
heavy  doping  (NHD),  bandgap  narrowing  (Slotbcom-DeGraaff 
model)  and  Fermi  Dirac  statistics  (BGN  +  FD),  bandgap 
narrowing  (Slotbocm-DeGraaff  model)  only  (BGN^,  and  for  a 
flat  profile  (NHD  +  flat) 


If 
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Sp  » 


Np(WE) 

n2  (W_) 
ie  E 


,Wr 


(5.21 


0 


Np(x) 
n2e(x) 


dx 


then  (5.20)  reduces  to 


TO 


Wr 


'Nq(x) 

n2  (x) 

ie  ' 


(5.22) 


dx 


which  is  the  exact  expression  for  an  infinite  surface  recombination 
velocity.  For 


Sn  « 


N0(WE) 


p    ZK) ' 


ie  t 


(5.23) 


n?e(x) 


dx 


(5.20)  reduces  to 


qSpnfe(WE) 


UP0 


N0rwE, 


[5.24) 


It  is  desirable  to  make  Jp_  small  for  the  bipolar  transistor. 
This  results  in  a  large  emitter  efficiency.  For  the  pn-junction  solar 


cell,  if  the  emitter  recombination  current  Jpn  is  small  compared  with 
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the  base  recombination  current,  the  value  of  V„f  can  approach  the 
classical  theoretical  limit. 

To  illustrate  the  dependence  of  Jp„  on  Sp,  consider  the 
desirable  case  in  the  transparent-emitter  model  in  which  Sp  is  small 
enough  to  satisfy  (5.23).  Figures  5.3  and  5.4  show  the  variation  of 
the  emitter  saturation  current  density,  JpQ,  and  the  transit  time,  x., 
as  a  function  of  Sp  for  three  models  of  bandgap  narrowing:  Slotboom 
and  De  Graaff  [20],  Hauser  [12],  and  the  recent  model  of  Lanyon  and 
Tuft  [13],  which  has  the  form 

AEG  =  22.5  x  10"3  (N/1018)1/2  eV 

for  non-degenerately  doped  silicon,  and  (5.25) 

AEQ  =  162.  x  10-3  (N/1020)1/6  eV 

for  degenerately  doped  silicon.  A  Gaussian  impurity  profile  is  assumed 
with  a  surface  impurtiy  concentration  of  1020  cm"3  and  a  junction  depth 
of  .25  ym.  Full  ionization  of  the  impurity  atoms  is  also  assumed. 

For  values  of  Sp  below  10  cm/sec,  Jpq  and  i\  vary  rapidly  with 
variations  in  Sp,  while  for  values  of  Sp  above  10°  cm/sec,  both  Jpp 
and  tt  saturate.  JPq  saturates  to  its  largest  value,  and  xt  saturates 
to  its  lowest  value.  The  largest  value  of  Jp0  at  any  Sp  occurs  for 
the  Lanycn-Tuft  model  of  bandgap  narrowing. 

The  validity  of  the  transparent-emitter  model  is  based  on  the 

condition  that  the  minority-carrier  transit  time  is  much  smaller  than 

the  minority-carrier  lifetime:   x+.  <<  t   .     To  test  this  condition,  x. 

t  p  t. 


53 


— 
o 

+-> 

n3 

c 


3< 


(UT3J3    O 

S-    (U    ra  +J 

■D 

<JJ    3    E    0) 

Ur-     fl     3 

re    (J  JS  "O 

4-    C  +-> 

S-  -i-           C 

3           on    0) 

n   ouir 

E    <1)   -l-> 

OJ   •!-   r— 

.e   3  ■---•  oi 

o 

S-     D. 

o 

4-  i_  O0  +J 
Ore           C 

th 

C  4-    0) 

c        o  s- 

*s 

O    Q.         S- 

£ 

■r  fl  1/1  3 
O  "O    3 

(  ) 

C    C  r-    S- 

3  (0  re  <u 

4—  .a    >  +-> 

(O-O    2-t- 


C/} 


oi   s_ 

a  > 

ou  o 

-c  o 

Q_-i-  U_ 

+->  -0 

■-> 

+-> 

re 

1/1 

o 

>! 

ul  TJ 

+J 

■+->  4- 

C1J 

n3  4— 

Q-,-4-> 

00 

+->    n3 

H    S- 

C 

i/l    ra 

O 

a) 

s_ 

a   a. 

"D 

U  CD 

0) 

re  oj 

•MS- 

j_> 

S_  Q 

H 

c 

<s> 

O) 

a  e 

1/1     OJ 

S- 

1    o 

■a  3 

s- 

•r-    O 

3 

=  -Q 

aj  re 

c J 

i.  -M 

•r-    > 

0)   o 

>> 

c 

u_ . — 

C 

o 

oo 

+->  re 

on  .c 

4-> 

^o 

0)   +J 

ra 

+-> 

<_ 

U-) 

S- 

3 

c\j 

>>  oj 

re  o  c  c  s- 


r- 

£. 

-s 

o 

> 

'3 

O 

iT) 

u 

re 

OJ 

4- 

z 

i 

jr 

x:    q-cq  a>  c 


(_uuo/v)odr 


59 


C/) 


.-    o 


c 

C 

o 

•  i— 

•  1 — 

z 

+-> 

o 

u 

s- 

c 

s_ 

3 

(0 

4- 

c 

« 

o. 

m 

l/> 

en 

(O  X3 

(03S)  »a 


50 


is  calculated  from  (5.11).  Values  of  Tx  are  plotted  in  Figures  5.4 
and  5.5.  In  Figure  5.4  xt  is  plotted  as  a  function  of  Sp  for  the 
three  models  of  bandgap  narrowing  (assuming  WF  =  .25  pm),  while  in 
Figure  5.5,  values  of  x.  are  plotted  as  a  function  of  the  quasi-neutral 
emitter  region  width,  W^,  for  Sp  =  5  x  10^  cm/sec.  The  recombination 
lifetime  x  has  an  upper  bound  determined  by  the  Auger  band-to-band 
recombination  at  high  impurity  concentration. 

To  illustrate  the  self-consistency  test  for  transparency,  we 

20    3 
assume  the  surface  concentration  of  a  diffused  emitter  to  be  10   cm  . 

With  the  impurity  profile  assumed  Gaussian,  this  corresponds  to  an 

average  Auger  lifetime1  of  x.  =  2.4  x  10   sec.  In  Figures  5.4  and 

5.5,  we  compare  this  lifetime  with  x,  for  each  of  the  three  bandgap 

narrowing  models  (assuming  that  WF  =  0.25  ym) .  Note  that  the  emitter 

5 
is  completely  transparent  if  Sp  exceeds  10  cm/sec  and  is  opaque  if 

4 
Sp  is  below  10  cm/sec. 

5.5  Application  to  pn-Junction  Silicon  Solar  Cells 
In  this  section,  we  apply  the  transparent-emitter  model  to 

calculate  the  open-circuit  voltage  of  silicon  pn-junction  solar  cells 

having  low  substrate  resistivity. 

The  open-circuit  voltage  is  given  by 

v   ,  kI£n^SC  (5>26) 

0C   q     JQ 


'The  average  Auger  recombination  lifetime  was  calculated  from  a 
model  similar  to  that  of  W.  W.  Sheng  [37],  in  which  the  emitter  is 
divided  into  two  regions;  one  region  has  Shockley-Read-Hall  recombina- 
tion, and  the  surface  region  is  dominated  by  Auger  recombination. 
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Figure  5.5  T+   versus  UV  for  S  =  5  x  10  cm/sec.  Fermi-Dirac 
statistics  and  bandgap  narrowing  are  included: 
(A)  Lanyon-Tuft,  (B)  Hauser,  and  (C)  Slotbocn-DeGraaff 
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where  J<-p  is  the  short-circuit  current  density  and  Jg  is  the  saturation 
current  of  the  solar  cell  in  the  dark.  The  saturation  current  density 
Jq  of  the  diode  is 

J0  =  °P0  +  °N0  (5'27; 

where  J   is  the  emitter  minority -carrier  saturation  current  density 
and  Jjvjq  is  the  base  minority-carrier  saturation  current  density. 
The  base  saturation  current  density  is 


2 

JMn  =  ir\±  (5.28 

NO   N^LN 


Consider  low-resistivity  silicon  solar  cells  with  base  doping  concen- 
tration of  Nnn  =  5  x  10  '  cm  .  Measurements  made  on  such  cells 
indicated  the  minority  carrier  diffusion  length,  L,,  to  be  80  urn  [66], 
corresponding  to  J»,q  -  6.2x10"'^  A/cm  .  The  general  expression  for 
the  transparent-emitter  current  density,  given  in  (5.8),  has  yielded 
Jpo  as  a  function  of  Sp,  as  shown  in  Figure  5.3. 

Combining  these  characterizations  for  Jmq  and  Jpg  with  J__  = 

'oc 


23  mA/cm  (AMO  conditions)  [1],  we  plot,  in  Figure  5.6,  VQC  versus  Sp. 
For  low  Sp  (about  10  cm/sec),  V   is  limited  by  the  base  current,  and 

for  higher  values  of  Sp,  VQC  is  limited  by  the  emitter  current,  as 

5 
has  been  observed  experimentally  [IJ.  Note  that  for  Sp  >  10  cm/sec, 

VQC  saturates  to  its  lowest  value.  Note  also  that  the  Lanyon  and 

Tuft  model  of  bandgap  narrowing  gives  lower  values  of  VgC  (for  any 

given  value  of  Sp)  than  those  given  by  the  Slotboom  and  De  Graaff  and 

Hauser  models. 
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So  far  we  have  assumed,  for  simplicity,  that  all  of  the  emitter 
surface  is  characterized  by  a  single  value  of  Sp.  He  now  consider 
a  more  realistic  structure  of  silicon  solar  cells,  Figure  5.7.  The 
emitter  saturation  current  Jpn  is  the  sum  of  three  components  from 
regions  1,  2,  and  3  shown  in  Figure  5.7(b),  [67].  The  components  of 
the  current  denisty  from  the  metal-covered  surface,  region  1,  and 
the  nonmetal-covered  surface,  region  2,  are  given  by  (5.8).  In  region 
3,  the  flow  of  minority  carriers  is  two-dimensional  since  the  minority 
carriers  within  about  a  diffusion  length  from  region  1  are  much  more 
influenced  by  the  high  value  of  Sp  of  region  1  than  they  are  by  the 
relatively  low  value  of  Sp  of  region  2.  To  avoid  the  complexity  of 
two-dimensional  analysis,  we  make  the  first-order  approximation  that 
the  component  of  Jpn  from  region  3  is  essentially  the  same  as  that 
from  region  1  (Jp,  -  Jp-, )  because  Sp  of  the  unmetallized  surface  can 
be  made  orders  of  magnitude  smaller  than  Sp  of  the  ohmic  contact. 

The  emitter  saturation  is  then 

IP0  ■  ^  +  A3)  Jpol  +  A2JpQ2  (5.29) 

where  A,  is  the  metallized  surface  area  and  (A~  +  A.)  is  the  unmetal- 
lized surface  area  and  J    and  Jpg2  are  the  corresponding  currents. 
The  area  A  is  approximately  equal  to 

A3  =  2nLQLp  (5.30) 

where  n  is  the  number  of  metal  grid  lines  (or  fingers),  Lq  is  the 
length  of  the  grid  lines  (see  Figure  5.7),  and  A„  =  A^  -  (A  +  A-j), 
AT  being  the  total  area  of  the  cell. 
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(a) 


Sp(mstal) 


Seff  «Sp(  metal) 


S  (non-me(ai) 


(b) 


Figure  5.7  (a)  The  structure  of  a  pn- junction  solar  cell 

(b)  The  three  components  of  the  emitter  current: 
Jp-i »  Jpp »  and  Jp3 
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As  a  numerical  example,  let  A  =  4  cm1-,  L  =  2  cm,  n  =  6, 

T  b 

Lp  =  1  urn,  and  assume  10%  metal  coverage.  Then  A-i  =  0.4  cm  , 

2  ? 

A2  =  3.59  cm  ,  and  A3  =  0.0024  cm.  In  this  case,  A3  is  negligible. 

and 


kT  ^<:r 

V0C  "  T  £n  A^ AT (5.31) 


2 
A^  *  JP01  +  ft".  *  JP02  +  JN0 


This  expression  can  be  used  to  estimate  Sp  of  the  nonmetall ized 
surface  from  experimental  values  of  Vq^.  For  diffused,  thin-junction 
pn-junction  solar  cells  made  on  low  resistivity  (  =  0.1  ft- cm)  material, 
the  maximum  observed  open-circuit  voltage  is  about  600  mV  [1].  As 
one  example,  if  we  consider  the  Lanyon-Tuft  model  of  bandgap  narrowing, 
and  let  Sp  of  the  ohmic  contact  be  10  cm/sec,  let  the  doping  concen- 
tration  be  Gaussian  with  a  surface  concentration  of  10   cm~J.  then  by 
using  (5.31),  with  A-j/Aj  =  0.1  and  VqC  =  600  mV,  we  get  S  (nonmetal)  = 
5  x  10  cm/sec.  Thus  the  value  for  S  could  result  in  the  low  Vqq 
seen  in  conventional,  diffused,  thin-junction  solar  ceils.  Note  that 
Figures  5.4  and  5.5  indicate  the  self-consistent  validity  of  the 
transparency  assumption  for  this  device  in  the  Lanyon-Tuft  model 
which  permits  use  of  (5.20). 

Although  the  preceding  discussion  has  focused  on  the 
transparent-emitter  model  applied  to  N  ?  silicon  solar  cells,  the 
model  can  also  be  applied  to  PTN  cells.  It  is  straightforward  to 
show  that  heavy-doping  effects  (bandgap  narrowing  and  Fermi-Dirac 
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statistics)  degrade  N+P  cell  performance  more  than  that  of  P  N  cells 
because  the  effective  mass  of  electrons  in  silicon  is  greater  than 
the  effective  mass  of  holes.  The  resulting  different  effective 
densities  of  states  in  the  conduction  and  valence  bands  (Nc  and  Nw) 
cause  the  onset  of  degeneracy  to  occur  at  lower  impurity  concentrations 
in  P-type  material  than  in  N-type  material  [32],  if  both  N-type  and 
P-type  regions  have  the  same  bandgap  narrowing.  Thus  the  net  effect 
of  bandgap  narrowing  and  Fermi -Dirac  statistics  is  to  degrade  the 
N-type  heavily  doped  region  more  than  the  P-type  region  with  the  same 
impurity  concentration.  This  may,  in  part,  be  responsible  for  the 
high  efficiency  P+NN+  cells  that  have  been  observed  experimentally 
[68-69]. 

5.6  Perspective 
This  chapter  has  dealt  with  the  transparent-emitter  model  of 
a  solar  cell,  which  is  defined  by  the  condition  that  the  minority 
carriers  in  the  dark  quasi-neutral  emitter  recombine  mainly  at  the 
surface  rather  than  in  the  bulk.  Surface  recombination  can  predominate 
over  bulk  recombination  if  the  emitter  junction  depth  is  shallow  enough 
and  if  the  surface  recombination  velocity  is  high  enough.  In  fact, 
this  occurs  in  typical  pn- junction  silicon  solar  cells,  as  demonstrated 
by  recent  experiments  showing  the  sensitivity  of  VqC  to  the  surface 
recombination  velocity  [29].  From  a  theoretical  standpoint,  the  self- 
consistency  test  in  Section  5.4  can  determine  the  validity,  for  a  given 
solar  cell,  of  the  transparent-emitter  model,  provided  the  emitter 
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recombination  center  density  is  low  enough  for  the  Auger  process  to 
dominate  over  the  Shock! ey-Read-Hall  process. 

Although  the  transparent-emitter  model  may  describe  many 
conventional  shallow  pn-junction  silicon  solar  cells,  the  high  value 
of  the  surface  recombination  velocity  S  necessary  to  validate  the 
transparent-emitter  model  is  not  necessarily  desirable  from  a  design 
point  of  view.  Growth  of  a  thermal  SiC^  layer  on  the  emitter  surface 
can  substantially  decrease  S  and  increase  Vqq.  For  such  devices,  the 
dark  emitter  recombination  current  is  determined  mainly  by  bulk 
recombination. 


CHAPTER  VI 

A  FIRST-ORDER  ENGINEERING-DESIGN  MODEL 
FOR  HEAVILY  DOPED  SILICON  DEVICES 


6.1  Introduction 

In  the  previous  chapter,  we  presented  a  rigorous  analytic  model 
for  the  emitter  that  included  bandgap  narrowing,  Fermi-Dirac  statistics, 
and  a  finite  surface  recombination  velocity  at  the  emitter  surface. 
The  main  assumption  of  that  treatment  is  the  transparency  of  the 
emitter  to  the  injected  minority  carriers,  that  is,  the  transit  time 
of  the  minority  carriers  is  much  less  than  their  average  bulk  lifetime. 
This  case  can  be  physically  realized  if  the  emitter  junction  depth  is 
shallow,  and  the  surface  recombination  velocity  is  high. 

In  this  chapter,  we  present  a  first-order  engineering-design 
model  that  includes  all  the  mechanisms  of  the  transparent  emitter  model 
except  that  we  relax,  in  this  treatment,  the  transparency  assumption 
and  include  a  position-dependent  Auger  lifetime.  Thus  we  provide  a 
general  first-order  model  that  includes  all  the  fundamental  heavy- 
doping  effects.  This  model  can  be  easily  used  to  characterize  the 
minority-carrier  current  in  the  heavily  doped  emitter. 

There  is  no  analytic  solution  for  the  minority-carrier  con- 
tinuity equation  when  a  position-dependent  lifetime,  impurity 
concentration,  and  n^e(x)  are  included.  To  provide  a  numerical  solu- 
tion, we  developed  a  computer  program  that  solves  the  continuity 
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equation.     This  computer  solution  serves  two  purposes.     First,   it 
provides  information  needed  to  characterize   the  first-ordar  model. 
Second,   it  verifies    the    accuracy  of  the  first-order  model. 

6.2  Modeling  Approach 

The  key  approximation  of  the  first-order  model  is  the  assump- 
tion of  a  negligible  net  drift  component  in  the  minority-carrier  current 
expression  in  the  heavily  doped  emitter  (see  Equation  (B.l)  in  Appendix 
B).  The  rationale  for  this  approximation  is  discussed  in  [29]. 

To  provide  further  justification,  we  note  that  the  approxi- 
mation is  consistent  with  calculations,  which  include  bandgap  narrowing 
and  Fermi-Dirac  statistics,  of  the  effective  doping  density,  i.e., 
(nf/nf  )N«g.  figures  6.1  and  6.2  snow  the  effective  doping  densities 
for  the  Lan>cn-Tuft  model  of  bandgap  narrowing  [13]  versus  the  acti'al 
dooing  densities  Nqq  and  N,,.  Note  that,  in  the  range  of  doping  from 
10   to  about  IC1  cm""1,  the  effective  doping  density  is  nearly 
independent  of  the  actual  doping  density,  which  corresponds  to  nearly 
zero  net  electric  field  for  minority  carriers. 

In  the  first-order  model,  therefore,  the  minority  earners  flow 
primarily  by  diffusion.  We  characterize  the  position-dependent  dif- 
fusion length  for  minority  carriers  by  an  average  value  L^: 

^  =  [D  •  ?A]1/2  (5.1) 

where  fT  is  an  average  value  for  the  Diffusion  coefficient  and  t  is 
an  average  Auger  lifetime  for  the  minority  carriers.  We  will  later 


II 


ro 


E 

o 


a 
a 


o 

u 

— 

•f— 

00  "3 

0) 

0)  -0 

a. 

3 

>i 

. — 

4-> 

O 

1 

c 

— 

••- 

i. 

0; 

O 

<♦- 

re 

aw 

GO 

'. -^ 

■r— 

--> 

^ 

'/l 

4J 

-, — 

— ' 

C/l 

ro 

- 

+-> 

V 

V> 

"O 

'.J 

-T-. 

m 

c 

'_ 

t/i 

V 

3 

O 

'/1 

s- 

c_ 

i. 

a; 

<T3 

> 

C 

tr-CL 

4~<a 

ojtr-. 

era 

QC 

^~ 

rc 

n 

>> 

1 

+j 

-> 

rr. 

(.1 

',- 

<- 

n> 

CU 

cz 

~^ 

0 

C.M- 

C 

0 

Ci. 

0 

0 

T3  "O 

0 

Q.' 

-=: 

> 

(£.tuo)^aaN 


73 


to  -o 

a 

Q  -O 

Q.  3 

>>>— 

■*->    O 

1     c 

Cl.  •<- 

s-  a 

o  s- 

1-    TO 

<Coi 

«a:u 

Z  -r- 

+-> 

>^   (/i 

+-:  t- 

•r-  -;-> 

o   <s 

c:  +■> 

a  oo 

T3 

<J 

cn  r3 

C    £.. 

o 

o-a 

1 

O     I 

TJ  -r- 

E 

i—    C 
(T3    <D 

o 

3  U_ 

U  T3 

to    C 

r3 

< 

<D 

< 

,M    = 

/i    2 

X 

3    O 

00     S_ 

S_    S_ 

0>    n3 

14—  Q_ 

<+-<T3 

(UP, 

<-3 

<c 

2:  ra 

.c 

>i   i 

+->   >^ 

—     CI 

0T    i- 

£=     O 

0)   c: 

-a  <u 

C7><4- 

c  o 

(s-UJ0)  "svvn 


4- 
fj 


74 


verify  that  the  use  of  L«  is  consistent  with  the  results  of  the  computer 
solution. 

6,3  Development  of  the  Mode"! 
Using  the  approximations  of  Section  6.2,  we  solve  the  minority- 
carrier  continuity  equation  in  the  heavily  doped  emitter  and  obtain 
the  following  expression  for  the  minority-carrier  saturation  current 
(Jpn  in  an  N+-emitter,  for  example): 


q  n7  Dp    SPLA  +  Dp  tanh(WE/LA) 


?0   NDDeff  '  CA   Dp  +  SpLA  tanh(WE/LA) 

where  Sp  is  the  emitter  surface  recombination  velocity  for  holes,  Wc 
is  the  emitter  quasi -neutral  region  width,  and  Ngrjeff  "is  an  average 
effective  doping  density. 

We  can  simplify  (6.2)  by  considering  the  following  special 
cases: 

A.  If  W-  »  La,  then  (6.2)  reduces  to 


(6.2) 


Unn 

r-U 


1  p 


"DDeff  '  LA 


(6.3) 


which  is  the  familiar  expression  of  the  minority-carrier  saturation 
current  for  a  completely  opaque  emitter. 

B.  If  F»  »  ME,  then  (6.2)  reduces  to 


TO 


q  nf    [1  +  WE/TASp] 
NDDeff   [l/3p  +  Wc/Dp] 


(5.4) 


We  further  simplify  the  expression  in  (6.4)  by  considering  different 
values  of  S_.  If  Sp  -  0,  then  (6.4)  reduces  to 


2 

q  n^  d£ 

Jon   -   = '  =~                                    (6.5] 

JDDeff  TA 


'P0   N| 


In  this  case  JpQ  is  directly  proportional  to  WE.  If  Sp  is  very  high 
(Sp  -*-™),  then  (6.4)  reduces  to 


q  ni    Dp 

jpc  -  j?    *  w7 

lMDCeff    L 


(5.5) 


and  Jpn  is  inversely  proportional  to  Wr.  For  the  case  when 


x„  »  tt^      and      Sn  << 

>P 


lA  "  Sp      """      ~P  ^  WE 

then  (6.4)  reduces  to 
2 

,]     =: L_  »  S  /'A  7' 

JP0      M  P  ^6"7. 

NDDeff 
and  JpQ  does  not  depend  on  W^. 

6.4     Relating  Device  Design  Parameters 
to  the  First-Order  Model's   Parameters 

In  this  section,  we  provide  means  for  characterizing  the 

parameters  of  the  first-order  model  for  any  set  of  design  parameters 

for  the  emitter.  At  least  three  design  parameters  need  to  be  known 

to  evaluate  the  recombination  current  of  the  emitter:  the  surface 
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impurity  concentration  N5,  the  emitter  quasi-neutral  region  width  UE, 
and  the  surface  recombination  velocity  S.  The  parameters  of  the 
first-order  model  are  Ne-Ff.  D,  tn,  WE,  and  S. 

We  now  describe  hew  PTef*»  D>  and  r,  are  determined.  We 
assume  that  the  impurity  profile  is  Gaussian.  The  parameter  NL^f  is 
obtained  from  a  spatial  average  using  Figures  6.1  and  6.2.  In 
Figures  5.3  and  6.4,  we  plot  average  effective  and  actual  doping 
densities,  NpQgff  and  N'DD  for  N-type  and  N^/uff  and  N.«  for  P-type, 
as  functions  of  the  surface  impurity  concentration  N$. 

To  find  D,  we  first  evaluate  the  average  doping  density  in 
the  emitter  (by  averaging  over  the  Gaussian  profile),  and  find  from 
taoles  the  mobility  of  the  minority  carriers  corresponding  to  N". 
By  using  the  Einstein  relation  in  it:  non-degenerate  or  degenerate 
[70]  form  (depending  en  N)  we  car.  calculate  D~. 

We  characterize  t^  as  a  function  of  N$  by  Figures  6.5  and 
6.6  for  N-type  and  P-type  silicon,  respectively.  These  characteriza- 
tions are  derived  from  the  rigorous  computer-aided  numerical  solution. 

6.5  Verification  of  First-Order  Mode' 


by  Computer  Analysis 
A  computer  program  has  been  developed  that  solves  the 
minority-carrier  continuity  equation  with  heavy-doping  mechanisms 
such  as  energy-bandgap  narrowing,  Fermi -Oirac  statistics,  and  a 
position-dependent  lifetime  included.  This  program  (see  Appendix  D 
for  a  listing)  is  used  to  verify  the  accuracy  of  our  first-order  model 
end  to  characterize  the  average  bulk  Auger  lifetime  in  the  emitter 
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Figure  6. 


The  average  actual  doping  density  Nqq  and  the  average  ef- 
fective doping  density  NpDeff  varsus  the  surface  impurity 
concentration  N$  for  N-type  silicon.  Lanyon-Tuft  model  of 
energy  bandgap  narrowing  and  Fermi-Qirac  statistics  a^e 
included. 
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Nc  (em3) 


Figure  6.4 


The  average  actual  doping  density  N^  and  the  average  ef- 
fective doping  density  NQDeff  versus  the  surface  impurity 
concentration  N5  for  P-type  silicon.  Lanycn-Tuft  model 
of  energy-bandgap  narrowing  and  Fermi-Dir-ac  statistics 
ere   included. 
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:igure  6.5  The  average  Auger  lifetime  r,  (calculated  by  computer- 
aided  analysis)  versus  N,.  for  N-type  silicon 
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Figure  6.6  The  average  Auger  lifetime  x^  (calculated  by  computer- 
aided  analysis)  versus  Nc  for  P-tyoe  silicon 
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as  a  function  o*   the  surface  impurity  concentration.  Sample  results 
of  the  computer-aided  analysis  aie  presented  and  discussed  in  Appendix 
E . 

In  Figures  5.7  and  6.8,  we  plot  the  minority-carrie1"  saturation 
currents  (and  their  surface  and  bulk  components),  J^q  in  a  PT-emitter 
and  JpQ  in  an  M  -emitter,  as  functions  of  SN  and  Sp,  respectively, 
for  Nc;  =  10"  cm   and  Wp  =  0.25  um.  i he  curves  are  derived  from  our 
computer-aided  analysis  and  the  points  are  derived  from  calculations 
using  the  first-order  mode"!  as  described  earlier  in  this  chapter. 
Note  the  good  agreement  between  the  results  of  the  first-order  model 
and  those  of  the  computer-aided  analysis. 

In  Tables  III,  IV.  and  V,  we  compare  values  of  the  minority- 
carrier  saturation  current  obtained  from  the  first-order  model  and  the 
numerical  analysis  for  N+-  and  P+-emitters  having  He  -   10^  cm"3 
(Table  III),  1020  cm"3  (Table  IV),  and  4  x  ID20  cm"3  (Table  V).  Note 
that  good  agreement  is  obtained  for  Nc;  =  1C   cm""  and  10L-L  cm"J  in 
both  N  -  and  P  -emitters.  For  Ms  =  4  x  10   cm  ,  fair  agreement  is 
obtained  for  the  N+-emitter,  but  poor  agreement  is  obtained  for  the 
?  -emitter.  The  inadequacy  of  the  first-order  model  for  this  latter 
case  can  be  explained  by  noting  the  differences  between  Figures  6.1 
and  6.2  In  Figure  6.1  (N  -emitter) ,  the  range  of  flatness  of  Mfjrjef  f 

on     -j 

extends  up  to  about  Ngg  =  4  x  i0dU  cm~J,  whereas  in  Figure  6.2 
(P  -emitter)  the  flatness  of   Nvicift-  extends  up  to  only  about  N*,  = 


L.      A 


10£U  erf  .  Hence,  at  about  N<-  -  4  x  10   cm  ,  the  first-order 


model  does  not  give  accurate  results  f'or  P  -regions,  although  it  gives 
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TABLE  III 

COMPARISON  OF  EMITTER  SATURATION  CURRENT 

CALCULATED  FROM  FIRST-ORDER  MODEL  AND  COMPUTER-AIDED 

ANALYSIS  FOR  Ns  =  10^9  cm"3-EXCELLENT  AGREEMENT 

IS  OBTAINED  FOR  BOTH  N+P  AND  P+N  DEVICES. 


WE=  0.25 

urn 

NS  = 

1019  cm-3 

N+P:    » 

IpgA/cm2 

P+N:   6 

NqA/oit 

S 

(cm/ sec) 

Numerical 
Solution 

First-Order 
Model 

Numerical 
Solution 

First-Order 
Model 

102 

1.9  x  10"14 

1.9  x  10"14 

8.2  x   10"15 

8.6  x  10"15 

103 

6.3  x  10"14 

6.7  x  10"14 

4.4  x   10-14 

5.0  x  10""14 

104 

4.5  x  10"13 

4.7  x   10"13 

3.8  x   10"13 

4.1   x  10"13 

105 

2.3  x  10"12 

2.3  x   10~12 

2.6  x   10"12 

2.9  x  10"12 

106 

3.3  x  10"12 

3.4  x  10"12 

6.1   x   10-12 

6.3  x  10'12 

10? 

4.1   x   10-12 

4.0  x   10-12 

7.0  x   10"12 

7.2  x   10"12 
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TABLE  IV 


COMPARISON  OF  EMITTER  SATURATION  CURRENT 

CALCULATED  FROM  FIRST-ORDER  MODEL  AND  COMPUTER-AIDED 

ANALYSIS  FOR  Ns  =  1020  cm'3— GOOD  AGREEMENT 

IS  OBTAINED  FOR  BOTH  N+P  AND  P+N  DEVICES. 


We  =  0.25  ym  Ns  =  1020  cm-3 

N+P:  JpoA/cm2  P+N:  JN0A/cm2 

S       Numerical     First  Order    Numerical     First-Order 
(cm/sec)     Solution       Model       Solution       Model 


102 

1.2  x  10"12 

9.9  x  10-'3 

1.2  x  10"13 

8.2  x  10"14 

103 

1.2  x  10-12 

1.0  x  10-'2 

1.3  x  10"13 

1.0  x  10"13 

104 

1.6  x  10"12 

1.4  x  10-12 

2.0  x  10"13 

3.1  x  10"13 

105 

3.2  x  10"12 

3.0  x  10"12 

8.0  x  10-13 

1.6  x  10"12 

106 

5.1  x  10"12 

4.5  x  10~12 

2.7  x  10"12 

3.4  x  10"12 

107 

5.4  x  1 0-1 2 

4.4  x  10"12 

3.8  x  10"12 

3.8  x  10-"12 

TABLE  V 

COMPARISON  OF  EMITTER  SATURATION  CURRENT 

FROM  FIRST-ORDER  MODEL  AND  COMPUTER-AIDED 

NUMERICAL  ANALYSIS  FOR  Ns  =  4  x  1020-F0R  THIS  CASE 

FAIR  ACCURACY  IS  OBTAINED  FOR  N+P  BUT 

POOR  ACCURACY  IS  OBSERVED  FOR  P+N. 


WE  =  0.25  urn  Ns  =  4  x  TO20  cm"3 


N+P:  JpoA/cm2  P+N:  JN0A/cm2 

S  Numerical  First-Order  Numerical     First-Order 
(cm/sec)     Solution       Model       Solution       Model 

102  1.9  x  10"12  6.5  x  10"13  6.9  x  10"14  1.9  xlO"16 

103  1.9  xlO"12  6.6  x  10"13  6.9x10-14  2.4  x  I0-16 

104  2.0  x  10"12  2.3  x  10"13  6.9  x  10"14  7.2  x  10"16 

105  2.0  x  10-12  i.o  x  10"12  6.8  x  10-14  3.9  x  10"15 
105  2.2  x  10"12  1.3  x  10"12  6.8  x  10"14  3.4  x  10"15 
107  2.3  x  10"12  1.3  x  10'12  6.8  x  10"14  9.6  x  10"15 


reasonably  accurate  results  for  N+-regions.  The  physical  reason 
for  the  differences  in  Figures  6.1  and  6.2  is  that  the  onset  of 
degeneracy  in  P+-emitters  occurs  at  a  lower  doping  density  than  for 
N  -emitters;  N^« e«  rises  sharply  due  to  degeneracy  of  the  charge- 
carriers  at  a  lower  doping  density  than  that  at  which  Nnr,  **   rises. 

6.6  Experimental  Support 
There  are  limited  experimental  data  that  support  the  first- 
order  model.  These  data  [29],  given  in  Table  II,  indicate  that  a 
reduction  in  the  emitter  surface  recombination  velocity  for  a  variety 
of  silicon  solar  cells  increases  the  open-circuit  voltage.  The 
emitter  currents  predicted  by  the  first-order  model,  and  verified 
by  the  computer-aided  analysis,  yield  values  of  Vgc  in  rough 
agreement  with  these  experimental  data.  We  will  consider  these  data 
further  in  the  next  chapter. 

6.7  Summary  and  Conclusions 
In  this  chapter,  we  established  a  first-order  engineering- 
design  model  for  heavily  doped  silicon  devices  and  solar  cells 
that  accounts  for  fundamental  heavy-doping  mechanisms.  After 
developing  the  model  and  discussing  its  assumptions,  we  related  the 
emitter  design  parameters  to  the  parameters  of  the  first-order  model. 
We  then  verified  the  accuracy  of  the  first-order  model  for  both  N+- 
and  P+-emitters  by  using  results  of  a  rigorous  computer-aided  analy- 
sis. We  found  that  the  first-order  model  is  valid  for  surface  con- 
centrations up  to  about  2  x  lO'-0  cm-3  for  P  -emitters  and  up  to  about 
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4  x  10^°  cm"  for  N+-emitters.  The  difference  in  the  range  of 
validity  of  the  first-order  model  between N  -  and  P+-  emitters  is 
attributed  to  the  onset  of  degeneracy  in  P+-  regions  occurring  at  a 
lower  impurity  concentration  than  in  N+-regions. 

The  first-order  model  is  useful  in  the  design  of  heavily 
doped  silicon  bipolar  devices  and  solar  cells.  The  model  provides 
a  simple,  but  sufficiently  accurate,  evaluation  of  the  emitter  re- 
combination current  that  can  be  effectively  used  in  such  designs. 


CHAPTER  VII 

COMPUTER-AIDED  STUDY  OF  V0c 
IN  N+P  AND  P+N  SILICON  SOLAR  CELLS 

7.1  Introduction 

In  this  chapter,  we  present  the  results  of  a  rigorous  computer- 
aided  numerical  analysis  of  the  heavily  doped  emitter  in  silicon  solar 
cells.  The  purpose  of  this  treatment  is  to  study  the  dependence  of 
Vqc  on  the  design  parameters  of  the  emitter  in  N+P  and  P+N  silicon 
solar  cells. 

We  include  the  following  fundamental  heavy-doping  mechanisms: 
bandgap  narrowing,  Fermi-Dirac  statistics,  and  Auger  recombination. 
We  also  include  a  finite  surface  recombination  velocity  at  the  emitter 
surface.  We  assume  low-level  injection  in  the  heavily  doped  emitter; 
this  assumption  does  not  significantly  limit  the  applicability  of  the 
analysis. 

7.2  Dependence  of  Vqq 
on  the  Emitter  Design  Parameters 
in  N+P  and  PTN  Silicon  Solar  Cells 


The  open-circuit  voltage  Vqq  is  given  by 


kT  ln    OSC 


'OC    q    JB0  +  JE0 


(7.1) 


where  Jsc  is  the  short-circuit  current  density  and  JgQ  and  JEq  are  the 
base  and  emitter  saturation  current  densities.  In  this  treatment,  both 


91 


92 


J5C  anc'  ^BO  are  determined  by  experimental  measurements  [1,66], 

and  Jrn  is  characterized  from  our  computer-aided  analysis  (see  Appendix 

E). 

In  studying  the  dependence  of  Vqq  on  J^q,   we  assume  that 
Jsc  and  Jgg  do  not  change.  This  is  a  reasonable  assumption  since 
variations  in  the  emitter  design  do  not  affect  Jgg  anc'  do  not  change 
Jsc  enough  to  cause  considerable  changes  in  Vqq.  We  let  J^q  = 
23  mA/cm2  [1]  and  JBQ  =  6.2  x  10"14  A/cm2  for  an  N+P  cell  and 
3.1  x  10~^4  A/cm2  for  a  P+N  cell,  which  correspond  to  a  diffusion 
length  in  the  base  of  about  30  ym  [66]  and  a  doping  density  of  about 
5  x  10"17  cm"3  for  N-type  and  P-type  substrates.  Note  that  our  assumed 
value  of  Jsc  is  not  essential  for  this  study  and  any  experimentally 
measured  value  can  be  incorporated  into  our  results  below  by  adjusting 
the  values  of  VQC  in  this  chapter  accordingly. 

7.2A  Dependence  of  Vqc  on  We 

We  have  studied  the  dependence  of  Vqc  on  W^  for  several  values 
of  the  surface  recombination  velocity  S,  in  both  N  -  and  P+-emitters. 
Two  values  of  the  surface  impurity  concentration,  N$  =  10   cm"°  and 
N$  =  1020  cm"  ,  have  been  used.  The  impurity  profile  in  all  cases 
considered  is  Gaussian. 

In  Figures  7.1  and  7.2,  we  plot  Vqc  versus  W^  for  values  of 
S  ranging  from  102  to  107  cm/sec  for  Ns  =  1019  cm"3.  These  curves, 
for  N+P  and  P+N  cells,  imply  the  following: 

A.  The  surface  recombination  velocity  significantly 
affects  Vqc  in  both  N+P  and  P+M  cells.  Thus, 
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for  low  surface  concentrations,  a  low  surface 
recombination  velocity  will  result  in  a  high  value 
of  Vqc- 

B.  The  range  of  the  surface  recombination  velocity 

that  significantly  affects  Vqq  is  102  <   S  <   10^  cm/ 

2 
sec.  Reducing  S  below  10  cm/sec  does  not  signifi- 
cantly affect  Vqq. 

C.  In  the  range  of  S  described  above,  S  affects  Voc 
more  in  P  N  cells  than  in  N  P  cells.  That  is,  the 
spread  in  values  of  Vqq  due  to  different  S,  for  any 
value  of  Wc,  is  larger  in  P+N  cells  than  in  N+P 

cells. 

5 

D.  Except  for  S  >  10  cm/sec,  Vqq  is  nearly  insensitive 

to  variations  in  W^.  For  S  >  10  cm/sec,  VqC 
increases  with  Wr. 

E.  Even  when  Wg  >  1  ym,  V™  is  affected  by  variations 
in  S  for  both  N+P  and  P+N  cells. 

The  main  conclusion  suggested  by  Figures  7.1  and  7.2  is  that 

19    ? 
for  low  emitter  doping  concentrations  (less  than  or  equal  to  10   cm-0), 

the  most  important  design  parameter  is  the  surface  recombination 

velocity.  3y  reducing  S,  VqC  can  be  increased. 

Of)      o 

The  dependence  of  Vqc  on  WE  for  NS  =  ^   cm   "is  snown  in 
Figures  7.3  and  7.4  for  N+P  and  P+N  cells,  respectively.  These  results 
imply: 
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A.  Variations  of  S  in  the  range  102  5  S  <  107  cm/sec 
affect  V0c>  for  both  N+P  and  P+N  cells,  but  less 

significantly  than  for  the  cells  having  Ns  = 

,  19   -? 
10   cm  d. 

B.  There  is  a  considerable  difference  in  the  effect  of 
S  in  N+-emitters  compared  to  P+-emitters.  In  N+P 
cells  S  affects  Vqq  considerably  if  Wf_  is  less  than 
0.5  ym,  while  in  P  N  cells,  S  affects  Vqq  appreciably 
if  We  is  less  than  1  ym. 

C.  The  spread  in  Vqq  due  to  variations  of  S  in  the 
range  102  <  S  f  10'  cm/sec,  at  low  values  of  Wr,  is 
larger  in  P+N  cells  than  in  N+P  cells.  For  example, 
if  W£  =  0.2  ym,  the  spread  in  Vqq  is  about  100  mV 
for  P  N  cells  compared  to  about  50  mV  for  N+P  cells. 

D.  For  both  N  P  and  P'N  cells,  Vqq  1s  sensitive  to 
variations  in  Wp  for  low  and  high  values  of  S,  but 

is  nearly  insensitive  to  variations  in  Wg  for  S  in 

4    5 

the  range  10-10  cm/sec. 

E.  For  large  values  of  WE(>  0.5  ym  in  N+P  cells  and 

1  ym  in  P  N  cells),  Vqq  saturates  at  about  630  mV 
for  P+N  cells  and  610  mV  for  N+P  cells. 
We  will  discuss  the  significance  and  the  physical  interpretation 
of  the  above  results  later  in  this  chapter. 
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7.2B  Dependence  of  Vqc  on  N$ 

We  have  studied  the  variation  of  VgC  as  a  function  of  Ns  for 
several  values  of  S  in  N+P  and  P+N  cells.  The  width  of  the  emitter 
quasi-neutral  region  for  these  studies  was  assumed  to  be  0.25  urn- 
We  plot  in  Figures  7.5  and  7.6  Vqc  versus  N<-  for  N+P  and  P+N  silicon 
solar  cells.  We  observe  that: 

A.  Surface  recombination  velocity  affects  Voc  tf  Ns 
is  less  than  4  x  1020  cm"3  for  N+P  and  P+N  cells. 
The  sensitivity  to  S  is  most  pronounced  at  the  lowest 
surface  doping  concentrations. 

B.  For  high  values  of  S  (>  104  cm/sec),  all  the  curves 
of  Vqq  versus  N^  have  minima  at  about  Ns  = 

'oc 


19    o  , 

5  x  10   cm"-3,  and  Vnr  increases  (slowly  for  N  P 


cells  and  rapidly  for  P  N  cells)  as  Ns  increases 


above  5  x  10   cm  . 


20  „m-3- 


C.  For  high  values  of  Ns  (>  4  x  10^u  cm"J),  V0c 

saturates,  independently  of  S,  at  about  680  mV  for 
P+N  cells  and  at  600  mV  for  N+P  cells. 


7.2C  Dependence  of  Vqc  on  S 

In  Figure  7.7,  for  both  N+P  and  P+N  structures,  we  compare 

°0      "3 

quantitatively  the  dependence  of  Vqq  on  S  for  Ne  =  1011  cm""  and 
W^  =  0.25  ym.  As  observed  before,  variation  in  S  affects  Vqc  wore 
in  P+N  cells  than  in  N+P  cells.  Note  that  Vqc  versus  S  curves  saturate 
at  lew  S  (<  10  cm/sec).  For  the  P+N  cell,  Vqc  saturates  at  about 
670  mV,  while  for  the  N+P  cell,  V0c  saturates  at  about  610  mV. 


Figure  7.5  The  open-circuit  voltage  Vqc  for  an  N+P  silicon  solar  cell 
versus  the  surface  impurity  concentration  N5  for  several 
values  of  the  surface  recombination  velocity  Sp 
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Figure  7.6  The  open-circuit  voltage  Vqq  for  a  P  N  silicon  solar  cell 
versus  the  surface  impurity  concentration  for  several 
values  of  the  surface  recombination  velocity  S»j 
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7.3  Discussion 
From  our  computer-aided  numerical  analysis,  we  deduce  two 
main  conclusions: 

1 .  Vqp  depends  strongly  on  S  for  moderate  surface 
impurity  concentrations  (i.e.,  for  N5  < 

4  x  102C  cm"3) 

2.  There  is  a  considerable  difference  in  the  values 
of  VQC  for  N+P  and  P  N  cells  having  identical 
structures. 

Because  the  above  results  are  derived  from  a  rigorous  computer- 
aided  analysis  that  incorporates  fundamental  heavy-doping  effects,  they 
are  inherent  in  the  device  structure  and  material  and  are  independent 
of  the  fabrication  process  for  the  cell.  Actual  devices  may  not  show 
the  theoretically  predicted  values  for  Vqc  if  process-induced  effects 
are  important.  For  example,  if  deep  traps,  induced  by  the  processing 
or  due  to  residual  impurities  like  gold,  dominate  the  recombination 
process,  then  our  predicted  dependence  of  Vqq  on  design  parameters 
will  not  be  consistent  with  experiment. 

Physically,  the  strong  dependence  of  Vqc  on  S  in  both  N+P 
and  P  N  cells  implies  that  the  emitter  in  these  devices  is  transparent 
to  the  injected  minority  carriers.  In  that  case,  the  analytical  study 
of  Chapter  V  is  applicable.  From  our  computer-aided  analysis,  we 
verified  that  the  emitter  is  transparent  for  low  surface  impurity 
concentrations  (<  4  x  1020  cm"3)  and  high  S  (>  104  cm/sec).  For 
Nj  -  10'y  cm,  the  emitter  is  transparent  for  S  larger  than  103  cm/sec. 


no 


The  difference  in  the  theoretically  predicted  values  of  Vnr 
for  N  P  and  P  N  cells,  with  similar  structures,  is  due  to  three  factors: 

1.  The  Auger  lifetime  of  minority  carriers  is  larger 
in  P+-silicon  than  in  N+-silicon  having  the  same 
impurity  concentration.  That  is, 

(xA>P+  >  'taV  (7-2) 

because,  as  experimentally  determined  [36],  the 
Auger  coefficient  in  P-type  silicon  is  smaller  than 
that  in  N-type  silicon. 

2.  The  onset  of  degeneracy  of  the  charge  carriers  in 
P+-silicon  occurs  at  a  lower  impurity  concentration 
than  in  N  -silicon.  This  is  due  to  the  lower 
effective  density  of  states  in  the  valence  band 
than  in  the  conduction  band.  Degeneracy  of  the 
charge  carriers  tends  to  decrease  the  emitter  current, 
as  discussed  in  Chapter  V,  thus  increasing  Vqq. 

Note  that  Voc  in  P+N  cells  is  higher  than  that  in 
N+P  cells  of  similar  structure  (see  Figures  7.5 
and  7.6). 

3.  The  mobility  and  diffusivity  of  the  minority 
carriers  in  the  base  region  of  (silicon)  N+P  cells 

are  larger  than  those  in  P+N  cells.  Hence  the  base 

+  + 

current  of  N  P  cells  is  larger  than  that  of  P  N 

cells,  if  both  cells  have  the  same  substrate  doping 
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concentration  and  the  same  diffusion  length.  In 
our  calculations,  we  had 

(JB0)N+p  «  2  (JB0)p+N  (7.3) 

The  effects  of  bandgap-narrowing  and  degeneracy  of  the  majority 
charge-carriers  can  be  seen  in  Figures  7.5  and  7.6.  As  N$  increases 
above  1019  cm-3  (for  S  >  104  cm/sec),  Vqc  decreases  in  magnitude, 
reaching  a  minimum  at  Ns  -  5  x  1019  cm-3.  The  decrease  in  v"oc  with 
increasing  Ns  is  mainly  attributed  to  bandgap  narrov/ing  which  increases 
the  emitter  current.  Above  Ns  -  5  x  10'9  cm"3,  degeneracy  of  the 
majority  charge-carriers  effectively  compensates  for  bandgap  narrowing; 
thus  Vqq  increases  with  increasing  Ns-  By  compensating  the  effect  of 
bandgap  narrowing,  degeneracy  enhances  the  built-in  electric  field 
due  to  the  impurity  concentration  gradient.  This  forces  the  minority 
carriers  away  from  the  surface,  and,  as  a  consequence,  Vqc  becomes 
less  sensitive  to  S.  With  the  minority  carriers  located  mainly  in 
the  low-impurity-concentration  region  (compared  to  the  surface  region) 
next  to  the  emitter  space-charge  region,  the  effective  Auger  lifetime 
is  not  so  short,  and  the  emitter  current  is  decreased  further.  At 
still  higher  values  of  Ns  (about  4  x  1020  cm"3),  Vqc  saturates,  tnat 
is,  increases  in  WE  or  S  do  not  affect  Vqc-  In  this  case,  the  confine- 
ment of  the  minority  carriers  away  from  the  surface  is  complete. 

Referring  back  to  Figures  7.5  and  7.6,  it  is  important  to 
point  out  that  the  effect  of  degeneracy  may  be  exaggerated  in  our 
treatment.  The  reason  for  this  is  fundamentally  related  to  the  density 
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of  states  in  heavily  doped  semiconductors.  This  is  discussed  in  Appendix 
F  briefly.  A  full  treatment  of  this  is  beyond  the  scope  of  this  disser- 
tation and  is  suggested  for  future  research. 

7.4  Experimental  Support 
Limited  experimental  data  indicating  the  increase  in  VqC  in 
typical  silicon  solar  cells  as  S  is  decreased  [29]  are  consistent 
with  our  theoretical  results.  Furthermore,  our  analysis  indicates  that 
Vqq  in  silicon  solar  cells  can  be  in  the  range  of  values  observed 
experimentally  (typically  600  mV  or  less)  because  of  heavy  doping  in 
the  emitter. 

7.5  Conclusions 

In  this  chapter,  we  presented  results  of  a  parametric  study  of 
Vqp  in  silicon  N  P  and  P+N  solar  cells  based  on  a  rigorous  computer- 
aided  analysis  of  the  minority-carrier  transport  in  the  emitter.  The 
emitter  current  is  influenced  by  heavy-doping  mechanisms  that  conse- 
quently limit  Voc  in  silicon  solar  cells. 

There  are  two  main  design-oriented  conclusions  of  this  chapter. 
First,  it  is  beneficial  to  have  a  low  emitter  surface  impurity  con- 
centration only  if  S  is  low  (about  10^  to  10  cm/sec).  That  is, 
Vqq  can  be  increased  toward  its  limit  value  of  700  mV  by  moderately 

doping  the  emitter  and  making  it  thin  and  transparent  to  the  minority 

3 
carriers  if  the  surface  recombination  velocity  is  below  10  cm/sec. 

Second,  it  is  easier  to  suppress  the  emitter  current,  and  thus  optimize 

Vqq,  in  P  -emitters  than  in  N+-emitters.  This  second  conclusion  is 

in  agreement  with  experimental  observations  [69,71]. 
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The  physical  reasons  behind  our  second  conclusion  are  (a) 
the  Auger  lifetime  in  the  P+-emitter  is  longer  than  that  in  the 
N  -emitter,  (b)  the  onset  of  degeneracy  occurs  at  a  lower  impurity 
concentration  in  P-type  silicon  than  in  N-type  silicon,  and 
(c)  the  mobility  and  diffusivity  of  electrons  in  the  base  of  N+P  cells 
are  larger  than  the  mobility  and  diffusivity  of  holes  in  the  base  of 
P+N  cells,  meaning  that  the  dark  base  current  in  N+P  cells  is  larger 
than  that  in  P+N  cells. 

Limited  experimental  data  and  our  theoretical  results  are 
consistent.  This  consistency  supports  the  validity  and  accuracy  of 
the  first-order  model  described  in  Chapter  VI  since  the  computer- 
aided  analysis  is  in  agreement  with  the  first-order  model. 


CHAPTER  VIII 
SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

8.1  Summary 

In  this  dissertation,  we  investigated  the  device  physics  of 
heavily  doped  silicon  devices  and  used  the  results  of  this  investiga- 
tion to  establish  engineering  models  for  design. 

In  Chapter  I,  we  specified  the  origin  of  the  research  problem: 
the  excessive  recombination  current  observed  experimentally  in  heavily 
doped  regions  of  silicon  pn-junction  devices.  This  observation  is 
in  conflict  with  the  predictions  of  classical  pn-junction  theory. 
The  excessive  current  is  due  either  to  a  large  excess  minority-carrier 
charge  in  the  emitter  or  to  a  short  minority-carrier  lifetime  or  to 
a  combination  of  these  causes.  Various  heavy-doping  mechanism  were 
then  examined  to  determine  the  simplest  physical  model  for  heavily 
doped  regions  in  silicon  that  is  consistent  with  the  observed  excessive 
recombination  current. 

In  Chapter  II,  we  evaluated  a  rigorous  analytic  model  for 
a  diffused  emitter  that  included  Auger  recombination  only,  i.e.,  a  model 
that  neglected  other  heavy-doping  effects.  We  found  that  such  a  model 
cannot  give  the  low  values  of  VqC  in  silicon  solar  cells  and  of  the 
common-emitter  current  gain  in  silicon  bipolar  transistors.  Thus 
we  concluded  that  other  heavy-doping  mechanisms,  in  addition  to  Auger 
recombination,  must  be  included  in  the  analysis  of  these  devices. 
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In  Chapter  III,  we  discussed  recombination  via  a  position- 
dependent  defect  density.  We  pointed  out  that  such  a  recombination 
cannot  yield  low  values  of  Vqc  in  silicon  solar  cells  and  the  experi- 
mentally observed  sensitivity  of  Vq^  to  surface  treatment.  We  then 
concluded,  considering  also  the  results  of  Chapter  II,  that  the  exces- 
sive recombination  current  is  not  only  due  to  short  lifetimes,  but  also 
to  an  excess  minority-carrier  charge  storage.  The  excess  charge 
storage  is  due  to  an  excess  intrinsic  carrier  density  n^e  in  the  heavily 
doped  region  (or  deficit  impurity  concentration  Neff). 

In  Chapter  IV,  we  showed  that  the  magnitude  of  n^e  and  N^ 
needed  to  predict  the  excessive  recombination  current  cannot  be  due  to 
deionization  of  impurities  as  described  in  the  recent  models  of 
Heasell  [25]  and  Popovic  [45].  This  deionization  does  not  remove  the 
need  for  an  effective  bandgap  narrowing  to  obtain  agreement  between 
theory  and  experiment. 

In  Chapter  V,  we  presented  and  discused  the  transparent-emitter 
model  in  which  the  minority  carriers  injected  into  the  emitter 
recombine  mainly  at  the  surface  rather  than  in  the  bulk.  We  showed, 
via  a  self-consistency  test,  that  the  transparency  of  the  emitter 
prevails  despite  Auger  recombination  in  the  bulk  provided  that  the 
surface  recombination  velocity  is  high  (above  10     cm/sec),  WE  is  less 
than  0.5  ym,  and  Ns  is  less  than  4  x  1020  cm"  .  Furthermore,  we 
verified  that  for  an  N^P  silicon  solar  cell,  with  emitter  parameters 
consistent  with  the  above  values,  the  transparent-emitter  model  can 
explain  the  experimentally  observed  value  of  Vqq  of  about  600  mV. 
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In  Chapter  VI,  we  presented  a  first-order  emitter  model  that 
can  accurately  predict  the  emitter  recombination  current.  The  key 
assumption  involved  is  a  zero  net  effective  drift  field,  that  is,  the 
assumption  that  the  quasi-electric  field  cancels  the  built-in  field. 
The  first-order  model  was  verified  quantitatively  by  comparing  its 
results  to  results  of  a  rigorous  computer  solution  of  the  minority- 
carrier  transport  problem  in  a  heavily  doped  emitter.  Limited 
experimental  support  for  the  first-order  model  was  briefly  considered 
in  Chapters  VI  and  VII. 

A  parametric  study  of  VQC  in  heavily  doped  N+P  and  P+N  silicon 
solar  cells,  based  on  our  computer-aided  analysis,  was  also  presented 
in  Chapter  VII.  The  study  reveals  the  dependence  of  VQC  in  silicon 
solar  cells  (both  N+P  and  P+N)  on  the  emitter  design  parameters,  and 
led  to  conclusions  that  are  important  to  design. 

8.2  Accomplishments 
In  this  dissertation,  we  accomplished  the  following: 

A.  Identified  the  important  heavy-doping  mechanisms  and 
showed  that: 

1.  Auger  recombination  alone  cannot  explain  the 
experimentally  observed  values  of  Vqq  in 
silicon  solar  cells  and  the  common-emitter 
current  gain  in  silicon  bipolar  transistors 

2.  Inclusion  of  Fermi-Dirac  statistics  is  im- 
portant, and  it  tends  to  compensate  the 
effects  of  bandgap  narrowing 
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3.  Deionization  of  impurities,  as  described  by 
the  recent  models  of  Heasell  and  Popovic, 
is  not  an  important  mechanism  in  heavily 
doped  silicon 

B.  Modified  the  semiconductor  transport  equations  (in 
Appendix  A)  in  heavily  doped  silicon  incorporating 
asymmetry  of  bandgap  narrowing  and  degeneracy  of  charge- 
carriers.  We 

1.  Modified  Shockley's  auxiliary  relations  to 
include  asymmetry  and  degeneracy 

2.  Showed  that  n^e  depends  on  the  total  bandgap 
narrowing,  and  on  the  asymmetry  in  bandgap 
narrowing 

3.  Developed  an  expression  for  the  minority- 
carrier  current  in  a  heavily  doped  region 
that  accounts  for  asymmetry  in  bandgap 
narrowing  and  Fermi -Oirac  statistics 

4.  Presented  a  new  approximation  of  Fermi-Dirac 
integrals  of  order  1/2  (in  Appendix  C) 

C.  Established  a  first-order  engineering  model  for  heavily 
doped  regions  in  silicon  pn-junction  devices 
(particularly  solar  cells)  and  verified  its  accuracy 

D.  Established  an  efficient  computer  program  that  solves 
the  minority-carrier  transport  problem  in  heavily 
doped  silicon  devices  with  the  fundamental  heavy- 
doping  mechanisms  included 
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E.  Studied  the  dependence  of  VQC  on  the  emitter  design 
parameters  in  N+P  and  ?~LN  silicon  solar  cells.  We 
concluded  that: 

1.  For  Ns  <  4  x  1020  cm"3,  surface  recombina- 
tion is  important 

2.  The  emitter  current  can  more  easily  be 
suppressed  in  P+N  than  in  N  ?  cells  but  is 
more  sensitive  to  S  in  P+N  cells 

3.  At  high  Ns  (above  about  4  x  1020  cm-3),  Vcc 
saturates,  that  is,  it  becomes  independent 
of  S  and  W£,  at  about  680  mV  for  P+N  cells 
and  at  about  600  mV  for  N+P  cells 

4.  Our  theoretical  results  are  consistent  with 
available  experimental  data 

8.3  Scope  and  Limitations 
The  discussion  of  the  physics  underlying  the  heavy-doping 
effect  has  been  limited  in  this  work  mainly  because  many  aspects  of 
heavy  doping  have  not  yet  been  fully  understood.  A  complete  theoretical 
study  of  the  physics  in  heavily  doped  semiconductors  is  beyond  the 
scope  of  this  dissertation.  We  considered  only  what  we  believe  are 
fundamental  first-order  heavy-doping  effects  and  neglected  or  roughly 
approximated  other  second-order  effects  (such  as  the  dependence  of  the 
density-of-states  effective  mass  on  the  doping  concentration). 

One  of  the  main  assumptions  is  the  rigid-band  approximation, 
that  is,  the  assumption  that  the  change  in  the  energy  bandgap  does  not 
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alter  the  shape  of  the  density-of-states  function  from  its  conventional 
square-root  dependence  on  energy.  Although  the  real  density  of 
states  in  heavily  doped  semiconductors  is  not  parabolic,  experimental 
measurements  of  the  differential  conductance  in  heavily  doped  N-type 
and  P-type  metal-silicon  tunnel  diodes  indicate  that  the  rigid-band 
approximation  is  good  for  P-type  silicon  and  is  fair  for  N-type  silicon 
[72].  Other  experimental  measurements  in  heavily  doped  silicon,  includ- 
ing the  electronic  specific  heat  [56]  and  nuclear  magnetic  resonance 
[57],  indicate  that  the  density  of  states  in  the  vicinity  of  the  Fermi 
level  is  well  approximated  by  a  parabolic  function. 
Other  approximations  made  in  our  study  are: 

1.  The  density-of-states  effective  mass  for  electrons 
and  holes  is  independent  of  impurity  concentration 

2.  Average  values  are  used  for  the  mobility  and  the 
diffusion  coefficient  in  the  transport  equations 

3.  Low-level  injection  prevails  in  the  heavily  doped 
emitter 

4.  The  phosphorus  or  boron  atoms  are  completely  ionized 
for  doping  levels  above  10'°  cm 

5.  One-dimensional  models  adequately  represent  the 
solar  cell 

6.  A  Gaussian  Profile  adequately  approximates  actual 
emitter  doping  profiles 

8.4  Recommendations  for  Future  Research 
Some  of  the  approximations  listed  above  deserve  further  con- 
sideration. In  addition,  we  have  several  other  recommendations  for 
future  research: 
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A.  Problems  pertaining  to  electronic  characterization 
in  heavily  doped  silicon  of  the  density  of  states, 
energy-bandgap  narrowing,  and  electron  affinity. 
These  can  be  explored  directly  by  experiments  in- 
volving measurements  of  differential  conductance 

in  heavily  doped  metal-silicon  tunnel  diodes,  magneto- 
absorption,  and  photoemission  in  heavily  doped  silicon. 
In  addition  extensive  electrical  measurements  of  the 
emitter  current,  as  in  [22],  for  N+P  and  P+N  diodes, 
with  the  use  of  our  theoretical  characterization  of 
the  emitter  current  can  provide  an  average  value  for 
n?  in  an  inhomogeneous  emitter. 

B.  A  new  model  to  calculate  deionization.  We  propose, 
for  future  investigations,  a  model  for  the  impurity 
energy  levels  in  heavily  doped  silicon  formed  by 
integrating  the  experimental  results  of  Pearson  and 
Bardeen  [47]  with  the  theoretical  results  of  Morgan 
[73]  pertaining  to  the  impurity  band.  In  the  pro- 
posed model,  the  impurity  density-of-states  function 
is  Gaussian.  Its  median  value,  which  is  the  value 
at  which  the  Gaussian  fucntion  has  its  peak,  is  a 
function  of  the  impurity  concentration  as  given  by  the 
experimental  model  of  Pearson  and  Bardeen  [47]  rather 
than  being  fixed  at  the  impurity  level.  We  anticipate 
that  such  a  model  may  agree  well  not  only  with  the 
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results  of  Pearson  and  Bardeen,  but  also 
with  other  experimental  results  discussed  in 
Chapter  IV. 

C.  Extension  of  our  analysis  to  high-low  junctions. 
Our  analysis  and  modeling  can  be  extended  to  study 
the  highly  doped  side  of  high-low  junctions  in 
back-surface-field  (BSF)  and  high-low-emitter  (HLE) 
silicon  solar  cells. 

D.  Fabrication  of  special  device  structures.  Experi- 
mental measurements  on  special  devices  that  are 
carefully  designed  and  fabricated  to  explore  par- 
ticular mechanisms  (for  example,  surface  recombination] 
together  with  our  theoretical  modeling  can  provide 
significant  understanding  of  the  physics  underlying 
silicon  device  operation  and  performance.  One  such 
design  that  deserves  study  is  a  pn-junction  with 

an  epitaxial  emitter  on  a  low-resistivity  substrate. 
The  emitter  is  fabricated  by  epitaxy  so  as  to  control 
the  shape  of  the  impurity  profile.  The  impurity 
profile  in  the  structure  is  a  "reversed"  profile;  that 
is,  it  has  its  highest  value  near  the  junction  and 
its  lowest  value  at  the  emitter  surface.  Such  a 
device  would  have  an  emitter  current  that  is  very 
sensitive  to  the  surface  recombination  velocity  since 
the  impurity  profile  aids  the  injected  minority 
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carriers  in  reaching  the  emitter  surface.  The  sur- 
face recombination  velocity  can  be  varied,  for 
example,  by  growing  thin  thermal  SiC^,  then  removing 
it  and  depositing  Si-jN*,  as  in  [29]. 
E.  Experimental  comparison  between  P+N  and  N+P  devices 
having  identical  structures  and  doping  levels.  Ex- 
periments to  provide  data  in  addition  to  the  limited 
experimental  data  presented  in  this  dissertation  are 
needed  to  further  corroborate  our  theoretical  conclu- 
sion that  P  -emitters  can  yield  lower  emitter  currents 


than  N+-emitters 
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APPENDIX  A 

EXTENSION  OF  SHOCKLEY'S  AUXILIARY  RELATIONS 
FOR  HEAVILY  DOPED  PN-JUNCTION  DEVICES 


Analysis  of  semiconductor  devices  is  most  often  based  on 
Shockley's  differential  equations  [11]:  the  current  density  equations, 

Jn  =  qDNVn  +  quNnE  (A.l) 

Jp  =  -qDpVp  +  qyppt  (A. 2) 

the  continuity  equations, 


|a-  V  •  jn  +  q(G  -  U)  (A. 3; 


|R-  -7  -  j+  q(G  -  U)  (A. 4) 


Poisson's  equation, 

7  •  t  =  72Ej/q  =  I   (p  -  n  +  NDD  -  NAA)  (A. 5) 

and  the  auxiliary  relations, 

n  =  n.  exp[EpN  "  El|  (A. 6) 

1       kT 


fEI  -  EFP]  ..   7] 

p  =  n.-  expl 1  lA./j 

1      kT   I 
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which  relate  the  carrier  concentrations  to  the  corresponding  quasi-Fermi 
levels.  The  symbols  used  above  have  their  usual  meanings  and  are  de- 
fined explicitly  in  [19]. 

The  auxiliary  relations  (A. 6)  and  (A. 7)  provide  a  useful  con- 
nection between  the  electron  and  hole  densities  internal  to  the  device 
and  the  voltages  at  the  device  terminals  [19].  This  connection  has  been 
exploited  in  device  analysis,  particularly  in  computer-aided  numerical 
analyses  [17,74-82]. 

For  devices  with  heavily  doped  regions  —  for  example,  the 
emitter  regions  in  bipolar  transistors  and  solar  cells—experiment  and 
theory  have  been  brought  into  agreement  by  accounting  for  the  effect 
of  energy-bandgap  narrowing  [29,83].  This  is  done  by  substituting  n? 
for  n^  in  the  np  product;  n^e  is  proportional  to  exp(AEg/kT)  [83],  in 
which  aEq  is  referred  to  as  the  bandgap  narrowing. 

Papers  that  have  accounted  for  bandgap  narrowing  in  computer- 
aided  numerical  analysis  [17,80-82]  have  not  defined  the  pre-exponential 
factors  that  appear  in  the  general  forms  of  the  auxiliary  relations 
valid  for  heavily  doped  semiconductors. 

In  this  appendix,  we  demonstrate  that  bandgap  narrowing  is  not 
properly  accounted  for  by  exchanging  n-  for  n^  in  (A. 6)  and  (A. 7),  as 
was  done  in  [81,82],  except  for  an  uncommon  special  case.  This  special 
case  is  that  of  symmetric  bandgap  narrowing,  for  which  the  shifts  in 
the  conduction-  and  valance-band  edges  are  equal  in  magnitude—which 
rarely  occurs  [84]--and  for  which  Boltzmann  statistics  describe  the 
distribution  of  the  majority  carriers—which  is  often  a  poor  approxima- 
tion. In  this  appendix,  we  also  derive  the  appropriate  pre-exponential 
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factors  in  the  generalizations  of  (A. 6)  and  (A. 7)  by  accounting  for 
asymmetric  bandgap  narrowing  and  Fermi-Dirac  statistics.  The 
resulting  modified  auxiliary  relations  are  shown  to  be  consistent  in 
form  with  the  expressions  for  minority-carrier  current  used  previously 
in  analytical  treatments  of  heavily  doped  regions  [16,83]  in  which  the 
drift  component  is  influenced  by  a  quasi-electric  field  resulting 
from  a  spatial  dependence  of  the  bandgap  narrowing.  Furthermore,  we 
show  that,  for  a  degenerately  doped  semiconductor,  n^  ,  the  square  of 
the  effective  intrinsic  carrier  concentration,  depends  on  the  asymmetry 
in  bandgap  narrowing.  This  dependence  of  n^  on  asymmetry  complicates 
the  interpretation  of  experimental  measurements  of  n?e  and  of  AEq  in 
degenerately  doped  semiconductors. 

Asymmetric  bandgap  narrowing  in  heavily  doped  regions  is  due  to 
several  physical  mechanisms:  impurity-band  widening  [73],  band  tailing 
[15,85],  electron-electron  interaction  [86],  and  screening  effects  [12, 
13].  These  mechanisms  produce  unequal  shifts  in  the  intrinsic 
conduction-  and  valence-band  edges  because  they  affect  electrons  and 
holes  differently.  As  an  illustrative  example,  consider  impurity- 
band  widening  in  a  heavily  doped  semiconductor  region  where  discrete 
energy  levels  in  the  bandgap  spread  into  a  band  of  impurity  levels.  If 
we  consider  an  uncompensated  semiconductor,  the  impuriry  band  is 
formed  in  the  vicinity  of  the  majority-carrier  band  edge.  The 
impurity  band  thus  effectively  shifts  the  intrinsic  majority- 
carrier  band  edge  and  leaves  the  intrinsic  minority-carrier  band  edge 
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essentially  unchanged.  So,  in  this  example,  bandgap  narrowing 
resulting  from  impurity-band  widening  is  due  primarily  to  a  shift  in 
only  one  of  the  band  edges. 

To  begin  the  derivation  of  the  modified  auxiliary  relations, 
consider  an  intrinsic  or  a  lightly  doped  semiconductor.  The 
intrinsic  Fermi  level  is  given  by 

Ec  +  Ev   !      Nv 

Ei=JV1  +  rkT^  (A-8) 

where  Ec  and  Ey  are  the  conduction  and  valence  band  edges  as  shown  in 
Figure  A-l,  and  Nc  and  Nv  are  the  electron  and  hole  effective  densities 
of  states.  For  this  case,  the  general  forms  of  the  auxiliary  relations 
simplify  to  (A. 6)  and  (A. 7),  in  which  EpN  and  Epp  are  the  electron  and 
hole  quasi -Fermi  levels. 

In  a  heavily  doped  semiconductor,  we  assume  that  the  intrinsic 
conduction  and  valence  band  edges,  Eq  and  Ey,  shift  to  EA  and  EA,  as 
illustrated  in  Figure  (A-l).  In  general,  we  can  relate  Eq   and  Ey  to 
Eq  and  Ey  as  follows: 

Ec  -  Ec  -  AAEQ  (A. 9) 

Ev  =  Ey  +  (1  -  A)AEQ  (A. 10) 

where  AEq  is  the  total  bandgap  narrowing,  comprised  of  a  component 
due  to  the  conduction-band  shift  AAEq,  and  a  component  due  to  the 
valence-band  shift  (1  -  A)aEq;  A  is  a  factor  between  zero  and  one,  whose 
value  depends  on  the  physical  mechanisms  underlying  the  bandgap 
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narrowing.  For  example,  in  an  N-type  heavily  doped  semiconductor,  if 
bandgap  narrowing  were  due  mainly  to  impurity-band  widening,  then  A  ~  1 

Consider  now  an  N-type  heavily  doped  semiconductor  with  an 
inhomogeneous  impurity  concentration;  a  similar  derivation  applies  for 
a  P-type  semiconductor.  Even  for  conditions  of  non-equilibrium,  we 
can  assume  that  the  electrons  (majority  carriers)  obey  Fermi -Dirac 
statistics,  while  the  holes  (minority  carriers)  follow,  to  a  good 
approximation,  Maxwell-Boltzmann  statistics,  provided  the  injection 
level  is  low  [11].  For  a  parabolic  density  of  states  in  the  conduction 
and  valence  bands  (rigid-band  approximation),  we  have 


n  =  NcF1/2 


fEFN  -  E'rl 


kT 


(A.11 


and 


P  =  Nv  exp 


fEV  "  EF?1 
kT   : 


(A. 12) 


In  (A. 11),  F-,  12   is  the  Fermi-Dirac  integral  of  order  1/2.  Combining 
(A. 9),  (A. 10,  (A. 11),  and  (A. 12),  we  obtain 


pn  =  NcNvF1/2 


EFN  -  Ec  +  AAEG)    fEv  +  (1  -  A)AEG  -  E^ 

exp  : 

{  kT     J   I        kT 


(A. 13) 


We  now  use  an  analytic  approximation  for  F-.  ,~(n)  (see  Appendix  C) 


F  M   -    exp(n) 
rl/2in;   l  +  C(n)  exp(n) 


:a.i4: 
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where  C(n)  is  a  function  of  n  described  in  Figure  (C-l).  Analytic 
approximations  for  C(n),  for  -4  <  n  <  +12,  which  introduces  errors  of 
less  than  2.5%,  are  given  in  Appendix  C.  Use  of  (A. 14)  in  (A. 13) 
yields 


pn 


2    fAEG 

ni  expbr 


exp 


EFN  "  EFP^ 


kT 


1  +  C(nc)  exp(nc) 


(a. is; 


where 


Hr  = 


EFN  -  EC  +  AAEG 
kT 


(A. 16' 


If  we  now  define  n.  by 
le  J 


2  fEFN  "  EFP| 


pn  =  n.  exp 

•e       kT 


(A. 17) 


we  get,  from  (A. 15)  and  (A. 17) 


2    (AEg 
n2  =   ni  exPlTT 


ie   1  +  C(nc)  exp(nc) 


(A. 18) 


for  an  N-type  semiconductor.  When  both  electrons  and  holes  follow 
Maxwell -Boltzmann  statistics,  n^  reduces  to  the  result  given  by  DeMan 
[33] 


ie   "i 


nn-  exp 


A£G 
kT 


(a. ig; 
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From  (A. 10)  and  (A. 12),  we  can  write 
(1  -  A)AEG' 


p  =  n-  exp 


kT 


exp 


EI  -  EFP! 
I   kT   j 


If  we  define  n-  by 


nip  =   ^  exp 


'(1  -  A)AEG| 


I    kT    J 


(A. 20) 


(A. 21) 


we  then  have 


,  :I  "  FP 
p  =  n.  exp' 

1P    I   kT   J 


ET  - 


(A. 22 


where  n^  can  be  thought  of  as  an  effective  intrinsic  carrier  density 
for  holes.  From  (A. 17)  and  (A. 20),  we  get 


n---       kT 


'iP    I 


(A. 23; 


We  therefore  define  nn-  by 

2 


n- 


n   n 


]_e 


(A. 24) 


so  that 


n  =  nin  exp 


fEFN  -  Ell 
kT 


(A. 25) 


where  n.  can  be  thought  of  as  an  effective  intrinsic  carrier  density 
in  3  

for  electrons. 
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From  (A. 18),  (A. 21),  and  (A. 24),  notice  that  n.  f   n.  f     n?  . 
Equations  (A. 22)  and  (A. 25)  are  the  modified  Shockley  auxiliary  rela- 
tions for  heavily  doped  N-type  regions. 


If  the  electrons  and  holes  are  non-degenerate,  n.jn  and  niD 


reduce  to 


nin  =  ni  exp| 


fAAEG] 


kT 


(A. 261 


and 


ni  =  ni  exp 


f(l  -  A)AEG 


kT 


(A. 27! 


J 


Thus,  even  for  non-degenerate  carrier  concentrations,  asymmetric  bandgap 
narrowing  gives  different  pre-exponential  factors  in  the  auxiliary 
relations  (A. 22)  and  (A. 25).  This  is  contrary  to  the  conclusions  of 
[35].  Note  that  the  pre-exponential  factors  are  equal  only  if  both 
electrons  and  holes  are  non-degenerate  and  if  bandgap  narrowing  is 
symmetric  (A  =  1/2) .  (The  bandgap  narrowing  is  symmetric  when  the 
shift  in  the  conduction  band  edge  is  equal  to  the  shift  in  the  valence 
band  edge. ) 

By  differentiating  (A. 22)  and  using  (A. 23)  with  the  assumption 
that  (dEpj^/dx)  -  0  in  the  (N-type)  heavily  doped  region,  we  find  that 


p   dx     p  q  In     dx 
ie 


n   dx 


qD  dp_ 
M  p  dx 


(A. 28) 
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The  right-hand  side  of  (A. 28)  is  the  hold  current  density  in  a  heavily 

doped  N-type  region  [16,83],  including  an  effective  drift  component 

dependent  on  n|e(x)  and  n(x)  (see  Appendix  B).  Thus,  (A. 28)  is  con- 
sistent with 


dE 


FP 


Jp  =  V  ~dx~~  <A-29* 

which  is  the  conventional  Shockley  expression  for  the  hole  current  in 
terms  of  the  quasi-Fermi  level.  Note  that  the  hole  current  density  as 
described  in  (A. 28)  depends  on  n-  as  given  by  (A. 18),  and  not  on  the 
factors  n.jp  and  n-jn.  Thus  when  Maxwell -Boltzmann  statistics  are 
applicable  for  both  electrons  and  holes,  the  minority-carrier  current 
in  (A. 28)  depends  on  the  total  bandgap  narrowing  irrespective  of 
asymmetry,  as  noted  in  [35].  However,  if  the  electrons  or  holes  are 
degenerate,  the  minority-carrier  current,  via  its  dependence  on  n?  in 
(A. 28),  depends  not  only  on  the  total  bandgap  narrowing  but  also  on  the 
asymmetry  in  bandgap  narrowing,  contrary  to  the  conclusion  in  [35]. 

We  have  shown  that  degeneracy  of  electrons  or  holes  and 
energy-bandgap  narrowing  in  a  semiconductor  requires  that  n.  depend 
on  the  total  bandgap  narrowing  and  on  the  asymmetry  in  bandgap 

narrowing.  This  result  is  particularly  important  in  interpreting  ex- 

2 
perimental  measurements  of  n^  and  of  AEg  in  degenerately  doped 

semiconductors.  In  (A. 18),  the  term  (1  +  C(ri(0)>  where  nr,  is  given 
in  (A. 16)  for  an  N-type  semiconductor,  can  be  considerably  larger 
than  unity  for  a  degenerate  semiconductor.  If  this  term  is  neglected, 
the  magnitude  of  bandgap  narrowing  deduced  from  experimental  measure- 
ments will  be  smaller  than  the  true  shrinkage  (in  the  rigid-band 


134 


model)  of  the  bandgap.  Further  complication  in  interpreting  experi- 
mental measurements  cf  n?e  arises  from  the  dependence  of  nc>  in  (A. 16), 
on  asymmetry  in  bandgap  narrowing.  Therefore,  detailed  interpretation 
of  experimental  measurements  of  n?e  anc|  0f  AEG  in  degenerately  doped 
semiconductors  requires  a  knowledge  of  the  asymmetry  in  energy-bandgap 
narrowing. 


APPENDIX  B 
QUASI-ELECTRIC  FIELDS  IN  HEAVILY  DOPED  SEMICONDUCTORS 

In  this  appendix,  we  show  that  the  minority-carrier  (hole) 
current  in  a  heavily  doped,  quasi-neutral  N-type  region,  for  low-level 
injection,  is  given  by 


Jn  =  q^nP 


fkT 

r  i 

q 

2 
> 

dn 


ie   1   dn 


dx    n   dx 


qD 


d£ 

p  dx 


:b.u 


We  also  show  that  the  effective  electric  field  acting  on  the  minority 
carriers  is  composed  of  two  components—one  related  to  the  position 
dependence  of  n?  ,  and  a  second  related  to  the  majority-carrier 
concentration  gradient. 

In  thermal  equilibrium,  the  principle  of  detailed  balance  re- 
quires that 


dP0 
jp  =  qupPgEp  -  qDp  -^-  =  0 


(B.2! 


and 


dn0 
Jn  =  qiiNn0EN  +  qDN  -^  =  0 


(B.3) 


where  Ep  and  E^  are  phenomenological  effective  electric  fields  acting 
on  holes  and  electrons.  These  effective  electric  fields  can  be 
expressed,  via  (B.2)  and  (B.3),  as 
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E  -  -E.  .  JL 

P     U    '  Pn 


and 


EN  = 


dPQ 
dx 


dnr 


(B.4) 


'0 


dx 


(B.5) 


For  non-equilibrium,  but  low-level-injection  conditions,  if  we 
assume  full  impurity  ionization  in  a  heavily  doped,  quasi-neutral  N- 
type  region,  then 


n  =  nQ  *  NDD 

Using  (3.6)  and  the  thermal-equilibrium  condition 

2 


p0n0  =  n.e(x) 
in  (B.4)  and  (B.5),  we  obtain 


(B.6) 


(B.7) 


EP=I 


and 


P  n 


J_ 

2 

ie 


dnie 
dx 


dN 


DD 


"DD 


dx 


:b.8] 


Em  = 


_N   _L 
^N  '  NDD 


dN 


DD 


dx 


(B.9) 


Note  that  the  expressions  for  Ep  and  E..  given  in  (B.3)  and  (B.9)  hold 
for  low-injection  as  well  as  for  equilibrium  conditions.  In  (B.8)  and 
(B.9),  we  see  that  En  is  not  equal  to  Ew.  Mote  that  E^  depends  on 
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the  space  variation  of  N  ,  while  Ep  depends  on  the  space  variation  of 

2  ? 

Nqq  and  n.  .  The  position  dependence  of  n.  is  due  to  the  position 

dependence  of  AE~  and  of  nc  as  given  in  (A. 18). 

The  relation  between  Ep  and  E  is  found  by  combining  (B.8) 

and  (B.9).  This  yields 


2 


1 


Dp   i    dn-    Mn 

i_ie  j 

Note  that  for  N  material  Dj^/yN  >  Dp/iip  =  kT/q  because  of  the  degeneracy 
of  electrons  [70,87]. 

If  we  substitute  (B.8)  into  the  non-equilibrium  form  of  (B.2) 
obtained  by  using  (3.6),  we  get  (B.l). 


APPENDIX  C 

INCLUSION  OF  DEGENERACY  IN  THE  ANALYSIS 
OF  HEAVILY  DOPED  SEMICONDUCTOR  DEVICES 


Fermi-Dirac  integrals  of  order  1/2  defined  by  [26] 

E^2dE  ... 

1  +  exp(E  -  ti)  [LA 


Fl/2(n)  ~  ji 


i    r      E1/2dE 


0' 


appear  in  the  theory  of  metals  and  semiconductors  because  the  charge 
carriers  obey  Fermi-Dirac  statistics  and  because  an  E1//2  dependence 
is  used  to  describe  their  density  of  states.  For  metals  and  semi- 
conductors, values  of  F-| /odl)  are  used  to  calculate  the  Fermi  energy 
(=  kTn)  from  the  concentrations  of  electrons  or  holes.  Often  Maxwell- 
Boltzmann  statistics,  which  reduce  (C.l)  to 

F1/2(n)  «  exp(n)  (C.2) 

are  used  for  semiconductors  because  electron  and  hole  concentrations 
are  not  as  high  as  in  metals.  However,  for  heavily  doped  semiconductor 
regions,  as  in  the  emitter  of  bipolar  transistors  and  solar  cells, 
(C.2)  introduces  a  large  error  and  (C.l)  must  be  used. 

The  Fermi-Dirac  integral  in  (C.l)  only  can  be  evaluated 
numerically;  values  of  F-j  /?(n)  are  given  in  Tables  [26,88,39]  for  n  typ- 
ically in  the  range  of  -4  to  +10.  Several  analytic  approximations 
for  F-|  ,0{r\)   are  found  in  the  literature  [24,26,90-96].  Most  of  these 
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approximations  require  short  intervals  of  n  over  which  different 
expressions  are  used  to  approximate  Fy2(n)-  This  makes  the  evalua- 
tion of  F-j ^(n)  tedious  and  inconvenient  in  computer-aided  analysis 
and  in  analytic  treatments  of  heavily  doped  semiconductor  devices  and 
possibly  also  in  ether  problems  in  physics. 

Recent  approximations  of  F1/2(n)  [95,96]  yielded  simple  and 
fairly  accurate  expressions  for  the  reduced  Fermi  energy  n  in  terms 
of  the  charge-carrier  concentration.  However,  these  approximations 
do  not  give  simple  expressions  for  F-|/2(n)  itself,  which  is  used  not 
only  to  calculate  the  Fermi  energy,  but  also  in  the  expression 
for  the  square  of  the  effective  intrinsic  carrier  concentration  n? 
in  a  heavily  doped  semiconductor.  The  square  of  the  effective 
intrinsic  carrier  concentration,  which  enters  into  most  of  the  ex- 
pressions of  the  electrical  parameters  in  heavily  doped  regions,  is 
defined  by 

n?e  '  Vo  =  Wl/2(T1C)F1/2<V  <C-3> 

where  Ng  and  Pq  are  the  equilibrium  concentrations  of  electrons  and 
holes  in  the  semiconductor,  Nc  and  Ny  are  the  effective  density  of 
states  for  electrons  and  holes,  and 

EF  "  EC 

nc=-V  (C-4) 

and 


-V  -  LF  ,   x 

nV  =   "I  (C.5) 

kT 
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In  this  appendix  we  present  a  simple  analytic  empirical 
approximation  for  the  Fermi-Dirac  integrals  of  order  1/2,  and  use  this 
approximation  to  derive  a  simple  analytic  expression  for  the  effective 
intrinsic  carrier  concentration  rn  in  a  heavily  doped  semiconductor 
having  energy-bandgap  narrowing  and  degenerate  majority-carrier 
distribution. 

As  in  [91],  we  express  F-| ^(n)  by 

sx  d  ( n ) 
F1/2(n)  -  y  ;  C(n)  exp(n)  (c-6) 

where  C(n)  is  a  function  of  n  given  graphically  in  [91].  In  Figure  C-l, 
we  have  calculated  C(ri)  from  tabulated  values  of  F-j/2(n)  [26].  Our 
Figure  C-l  is  similar  to  Figure  4  in  [91]  except  that  we  have  extended 
the  range  of  n  from  -4  to  an  upper  bound  of  +12. 

By  curve-fitting  techniques,  we  found  the  following  analytic 
expressions  for  C(n): 

C(n)  =  -4.4  x  10"2n  +  3.1  x  10"1    for  -4  <  n  <   +2       (C.7) 

and 

C(n)  =  exp(3.6  x  10"3n2  -  3.2  x  10_1n  -  8.8  x  10"1) 

for  +2  <  t)   _<+12      (C.8) 

These  expressions  for  C(n)  provide  values  of  F\/2   over  tne  range 

-4  5  n  -   +12  with  less  than  2.5%  error.  We  have  reported  an  expression 

similar  to  (C.4)  in  a  previous  paper  [24].  In  this  note  the  accuracy 
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of  (C.7)  is  increased  over  the  counterpart  in  [24]  so  that  the  error 
in  the  range  -4  <  n  <  +2  is  less  than  1%.  For  the  expression  in  (C.8), 
the  largest  error  lies  in  the  range  +2  <  n.  <  2.5  and  is  about  2.2%. 
For  values  of  n  >  2.5,  the  expression  in  (C.8)  yields  an  error  in 
F,  .  (n)  less  than  1%. 

In  Table  C-I,  we  present  values  of  F-|  nM   f°r  several  values 
of  n  in  the  range  -4  _<  n  5  +12.  Exact  values  of  F-i/odl)  are  taken  from 
tables  in  [25],  and  approximate  values  are  calculated  using  (C.6), 
(C.7),  and  (C.8).  The  percentage  error  is  shown  in  the  fourth  column 
of  Table  C-I.  Note  that  the  error  is  less  than  1%;  the  largest  error 
is  in  the  range  +2  <  n  <  +2.5  and  is  about  2.2%. 

If  we  use  (C.6)  in  (C.3),  we  get 

2 NCNV  exp(-EG/kT) 

nie  =  [1  +  C(nc)  exp(nc)]  •  [1  +  C(nv)  exp  (nv)]        (C'9' 

where  n_c  and  ny  are  given  in  (C.4)  and  (C.5),  and  Eq  is  the  energy- 
bandgap  of  the  heavily  doped  semiconductor, 

EG  =  EGQ  -  AEQ  (CIO) 

In  (CIO),  Eqq  is  the  energy-bandgap  of  the  intrinsic  (pure)  semi- 
conductor, and  AEp  is  the  bandgap  narrowing  due  to  high  impurity  con- 
centrations in  the  doped  semiconductor.  Equations  (C.9)  and  (CIO)  can 

2 
be  combined  with  the  definition  of  n.  [11]  to  give 

2  n^  exp(AEG/kT) 

Rie  =  [1  +  C(nc)  exp(nc)j  •  [1  +  C(nv)  exp(nv)]        ''"'" 
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In  (C.ll),  we  have  derived  a  simple  analytic  expression  for 

p 

n.  that  provides  physical  insight  into  the  effects  of  bandgap  narrow- 
ing and  degeneracy  of  the  charge  carriers  on  n?  .  Note  that,  for  the 
most  part,  the  effect  of  energy-bandgap  narrowing  is  included  in  the 
numerator  of  the  right-hand  side  of  (C.ll),  while  degeneracy  is 
incorporated  in  the  denominator  of  the  right-hand  side  of  (C.ll)  (for 
a  more  rigorous  analysis  of  the  dependence  of  n?e  en  degeneracy  and  bandgap 
narrowing,  see  Appendix  A).  From  (C.ll),  it  is  clear  that  energy- 
bandgap  narrowing  increases  n^e,  while  degeneracy  or   the  charge 
carriers  decreases  n^e  by  increasing  the  terms  in  the  denominator. 
If  only  one  type  of  carrier  is  degenerate,  as  is  usually 
the  case  for  the  majority  carriers  in  heavily  doped  semiconductor 
regions,  (C.ll)  reduces  to 

2     n?  exp(AEQ/kT) 
nie  =  1  +  C(nc)  exp(nc)  (CJ2) 

for  an  N-type  region,  and 


2     A   exP(AEG/kT)  .  u) 

nie  =  1  +  C(nv)  exp(nv)  [Lm     ' 

for  a  P-type  region.  Furthermore,  if  both  kinds  of  carriers  are  non- 
degenerate,  (C.ll)  reduces  to  [83]: 

n?  =  n?  exp(AEr/kT)  (C.14) 
ie    i      « 

In  conclusion,  we  presented  a  simple  analytic  approximation 

of  the  Fermi-Dirac  integral  of  order  1/2  that  is  suitable  for  analysis 
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of  heavily  doped  semiconductor  devices.  We  used  our  new  approximation 
to  derive  simple  analytic  expressions  for  n?  ,  the  square  of  the 
effective  intrinsic  carrier  concentration,  in  the  heavily  doped  semi- 
conductor. 
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TABLE  C-I 


COMPARISON  BETWEEN  EXACT 


AND  APPROXIMATE  VALUES  OF  F 


1/2 


In) 


n 

Exact 
Fn /?(n) 

[26] 

Approximate 
F1/?(n)  Using 
(6.6),  (C7), 

and  (C.8) 

Percentage 
Error 

-4 

1.82 

X 

io-2 

1.32  x 

IO"2 

0.2 

-3 

4.89 

x 

10-2 

4.88  x 

IO"2 

0.3 

-2 

1.29 

X 

IO"1 

1.29  x 

io-i 

0.5 

-1 

3.28 

X 

io"1 

3.27  x 

IO"1 

0.2 

0 

7.65 

X 

io"1 

7.69  x 

io-i 

0.5 

1 

1.58 

1.59 

0.1 

2 

2.82 

2.81 

0.9 

3 

4.49 

4.49 

0.1 

4 

6.51 

6.57 

0.9 

5 

8.84 

8.92 

0.8 

6 

1.14 

X 

10' 

1.15  x 

IO1 

0.6 

7 

1.43 

X 

IO1 

1.44  x 

10 

0.4 

8 

1.74 

X 

10' 

1.74  x 

IO1 

0.6 

9 

2.06 

X 

IO1 

2.08  x 

IO1 

0.7 

10 

2.41 

X 

IO1 

2.42  x 

IO1 

0.6 

11 

2.77 

X 

IO1 

2.78  x 

IO1 

0.2 

12 

3.15 

X 

IO1 

3.13  x 

IO1 

0.7 

APPENDIX  D 
LISTING  OF  THE  COMPUTER  PROGRAM 

The  computer  program  we  developed  for  our  analysis  is  listed 
below.  Parts  of  the  program  are  from  the  Harwell  Subroutine  Library, 
Harwell,  England. 
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APPENDIX  E 

CHARACTERIZATION  OF  THE  EMITTER  CURRENT 

IN  HEAVILY  DOPED  SILICON  DEVICES 

BY  COMPUTER-AIDED  NUMERICAL  ANALYSIS 


In  this  appendix,  we  present  characterizations  of  the  emitter 
current  as  a  function  of  design  parameters:  surface  impurity  concen- 
tration, surface  recombination  velocity,  and  the  width  of  the  quasi- 
neutral  emitter  region.  This  study  is  based  on  our  rigorous  computer- 
aided  analysis  of  the  minority-carrier  transport  in  a  heavily  doped 
semiconductor  region.  A  listing  cf  the  computer  program  is  given  in 
Appendix  D. 

In  Figures  E-l  and  E-2,  we  plot  the  variations  of  the  emitter 
saturation  currents  JP0  (for  an  N+-emitter)  and  J^q  (for  a  P+-emitter), 
including  both  the  surface  and  the  Auger  (bulk)  components,  as  functions 
of  S.  For  values  of  S  lower  than  104  cm/sec,  the  emitter  current  is 
mostly  due  to  the  Auger  (bulk)  component,  while  for  values  higher  than 
104  cm/sec,  the  emitter  current  is  almost  entirely  due  to  surface 
recombination. 

Figures  E-3  and  E-4  show  the  dependencies  on  S  of  xp  and  xM, 
the  effective  average  hole  and  electron  lifetimes  in  N  -  and  P  -emitters 
defined  by  the  charge-control  relation 
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where  Q£  is  the  total  charge  in  the  emitter,  and  JV  is  the  emitter 
current.  From  (E.l),  we  can  deduce  that 

i  =  ^  +  i  c«> 

where  xt  is  the  minority-carrier  transit  time,  and  xA  is  the  average 
Auger  (bulk)  lifetime  for  minority  carriers.  Note  that  for  S  larger 
than  10  cm/sec,  where  xt  <  xA,  the  emitter  is  transparent  to  the 
minority  carriers,  and  for  S  smaller  than  1CT  cm/sec,  the  emitter  is 
opaque. 

In  Figures  E-5  and  E-6,  we  plot  JpQ  and  J^q  versus  to-  for 
values  of  S  in  the  range  102-107  cm/sec  and  for  N-  =  1019  cm"3.  Note 
the  strong  influence  of  S  on  JpQ  and  J  -.  Note  also  the  weak  dependence 
of  JpQ  and  Jmq  on  Wjr.  For  this  case,  the  emitter  saturation  current 
density  can  be  approximately  given  by  the  first-order  model  current 
expression  in  (6.7) ,  in  which  the  emitter  current  is  directly  propor- 
tional to  S  and  is  independent  of  Wr. 

Figures  E-7   and  E-8  are  analogous  to  Figures  E-5  and  E-6, 

20    3 
but  with  N<-  =  10   cm-0.  Note  the  weaker  dependence  on  S  and  the 

stronger  dependence  on  W^  of  the  emitter  current  compared  to  the  vari- 
ations of  Figures  F-5  and  F-6.  Therefore,  we  conclude  that  as  N 
increases  above  10 'y  cm"  ,  the  emitter  current  becomes  less  dependent 
on  S  and  more  dependent  on  Wr,  provided  x+  remains  less  than  xa. 

Furthermore,  we  see  in  Figures  E-7  and  E-3  that  the  dependencies 
of  JpQ  and  J^q  on  W^arealtered  by  variations  in  S.  For  low  values 
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Figure  E-6  The  electron  saturation  current  density  JNq  in  a  P  -emitter 
versus  the  width  of  the  quasi -neutral  emitter  region  Wg  for 
several  values  of  the  electron  surface  recombination 
velocity  Sm 
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of  VJ£  (less  than  0.5  urn  for  N  -emitters  and  less  than  1.0  urn  for 
P  -emitters)  Jpg  and  J^q  vary  as  follows: 

A.  In  proportion  to  Wp_  for  low  S 

B.  Independent  of  WV  for  moderate  S,  and 

C.  In  inverse  proportion  to  Wp  for  high  S 

These  dependences  of  the  emitter  current  on  Wr,  for  different  values  of 
S,  are  consistent  with  the  expressions  of  the  emitter  current  given  by 
the  first-order  model  expressions  (6.5),  (6.6),  and  (6.7). 

In  this  appendix,  we  have  provided  plots  of  theoretical  results 

showing  the  variations  of  the  emitter  saturation  current  density  in 

+      + 
N  -  and  P  -emitters  as  the  emitter  design  parameters  vary.  These 

results  are  used  to  study  the  parametric  dependence  of  Vqc  "in  N+P 

and  P+N  silicon  solar  cells  on  design  parameters  in  Chapter  VII.  The 

results  in  this  appendix  can  also  be  used  to  characterize  the  emitter 

injection  efficiency  of  silicon  bipolar  transistors  once  the  base 

current  is  determined. 


APPENDIX  F 

ON  THE  PARABOLIC  DENSITY 
OF  STATES  IN  HEAVILY  DOPED  SILICON 


In  this  dissertation,  we  assume  the  rigid-band  approximation 
to  be  valid.  That  is,  we  assume  a  parabolic  density  of  states' 
dependence  on  energy,  i.e., 

g(E)=ceE1/2  (F.l) 

where  Cp  is  a  constant  coefficient,  independent  of  E.  In  this  appendix, 
we  discuss  the  applications  of  the  parabolic  density  of  states  to 
heavily  doped  semiconductors. 

High  impurity  concentrations  (typically  above  I019  cm"3  in  sil icon) 
introduce  a  large  number  of  "excess"  available  states  for  the  charge 
carriers.  Most  of  these  states  are  close  to  the  band  edge..  This 
increases  the  density  of  states  in  (F.l),  although  it  may  not  alter 
the  parabolic  nature.  In  other  words,  the  parabolic  density  of  states 
may  still  be  a  good  approximation,  as  suggested  from  experimental  ob- 
servations [56-57,72],  however  Ce  may  be  larger  than  its  value  in 
intrinsic  material.  If  Ce  increases  with  doping,  the  density  of  states 
remains  parabolic  but  the  concavity  of  the  parabola  describing  g(E) 
increases.  Hence,  there  are  more  available  states  than  in  the  pure 
semiconductor  at  any  energy  level  in  the  band.  This  mechanism  tends 
to  reduce  the  effect  of  carrier  degeneracy  as  considered  in  our 
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treatment.  We  cannot  comment  quantitatively  on  this  effect  without 
knowing  exactly  the  actual  density  of  states  in  heavily  doped  silicon 
and  fitting  them  into  parabolic  density  of  states  with  appropriate 
values  for  Ce. 

In  any  case,  degeneracy  of  the  majority  carriers  would  always 
have  a  lower  threshold  in  the  P+-emitter  compared  to  the  N+-emitter. 
Further  investigation  of  the  parabolic  approximation  for  the  actual 
density  of  states  is  suggested  for  future  research. 


REFERENCES 


[1]  F.  A.  Lindholm,  A.  Neugroschel,  C.  T.  Sah,  M.  P.  Godlewski, 
and  H.  W.  Brandhorst,  Jr.,  "A  Methodology  for  Experimentally 
Based  Determination  of  Gap  Shrinkage  and  Effective  Lifetimes  in 
the  Emitter  and  Base  of  p-n-Junction  Solar  Cells  and  Other  p-n- 
Junction  Devices,"  IEEE  Trans.  Electron  Devices,  ED-24,  402- 
410  (1977). 

[2]  A.  Neugroschel,  F.  A.  Lindholm,  and  C.  T.  Sah,  "A  Method  for 
Determining  the  Emitter  and  Base  Lifetimes  in  p-n  Junction 
Diodes,"  IEEE  Trans.  Electron  Devices,  ED-24,  662-671  (1977). 

[3]  H.  W.  Brandhorst,  Jr.,  Record  of  9th  Photovoltaic  Specialists 
Conf.,  1  (1972). 

[4]  W.  L.  Kauffman  and  A.  A.  Bergh,  "The  Temperature  Dependence  of 
Ideal  Gain  in  Double  Diffused  Silicon  Transistors,"  IEEE  Trans. 
Electron  Devices,  ED-15,  732-735  (1968). 

[5]  D.  Buhanan,  "Investigation  of  Current-Gain  Temperature  Dependence 
in  Silicon  Transistors,"  IEEE  Trans.  Electron  Devices,  ED-16, 
117-124  (1969).        "  ~ 

[6]  E.  0.  Johnson  and  A.  Rose,  "Simple  General  Analysis  of  Amplifier 
Devices  with  Emitter,  Control,  and  Collector  Functions," 
Proc.  IRE,  47,  407-418  (1959). 

[7]  R.  D.  Middlebrook,  "A  Modern  Approach  to  Semiconductor  and 

Vacuum  Device  Thoery,"  IEE  Proc,  106B,  Suppl .  17,  887-902  (1960). 

[3]  H.  K.  Gummel  and  H.  C.  Poon,  "An  Integral  Charge  Control  Model 
of  Bipolar  Transistors,"  Bell  System  Techn.  J. ,  49,  827-852 
(1969). 

[9]  F.  A.  Lindholm,  S.  S.  Li,  and  C.  T.  Sah,  "Fundamental  Limitations 
Imposed  by  High  Doping  on  the  Performance  of  PN  Junction  Silicon 
Solar  Cells,"  Record  of  the  11th  Photovoltaic  Specialists 
Conf.,  3-12,  1975. 

[10]  J.  G.  Fossum,  R.  D.  Nasby,  and  S.  C.  Pao,  "Physics  Underlying 
the  Performance  of  Back-Surface-Field  Solar  Cells,"  accepted 
for  publication  in  IEEE  Trans.  Electron  Devices,  April  1980. 

135 


136 


[12]     J.  R.   Hauser,   Final   Report  on  NSF  Grant  GK-1615   (1969). 

[13]     2o  SeaSily^d"?! -con-"TlEEc  :Band^P  Arrowing  in  Moderately 
1014-1078  (1975?  '    iSJ>ai^  Electron  Devices.  ED-26, 


[14]     Depende^fDenli-tv'o?  l±  '-indhol.n     "I.npunty  Concentration 
c.f:™:nJ  c  ??1^        States  and  Resulting  Fermi   Level   for 
Silicon,"  Solid, State  Electronics.   14,  4X7-416  (1971). 

[15]     V.   L.   Bonch-Bruyavich,  The  Electronic  Theory  of  Heavllv  iwh 
.Serni^c^idj^ors,  ^er^Tm^^FTT^Zl  York?  1966.     P~ 

[16]     R.  J     Van  Overstraeten,  H.  J.  OeMan,  and  R     p    Mertens     "Tr-r, 

K     qU3'n°onn  1n  Heavily  D°Ped  Si^on,"     EE     T^T  ElectroT 
Devices,   tD-20,   290-298  (1972).  d'i^_>iectron 

[1?]    theSdu*ent  Gl^T^nf1'0?5  <"  K«v"*  ^  Silicon  and 

['?]     5,vhVFair\"The  Effect  of  Strain-Induced  Band-Gap  Narrowina  on 

[19]     -.   A.   Lindholm  and  C.   T.   San,   "Fundamental   Electronic  Mechanic 
Limning  the  Performance  of  Solar  Cells,"   IEEE  Trans     Elect™ 
jevices,   ED-24,   299-304  (1977).  -  rans"   blect^on 

[20]     J     w.   Slotboom  and  H.   C.   DeGraaff,   "Measurement  of  Pandaan 
.arrowing  in  Silicon  Bipolar  Transistors.^olid  Stat Elec- 
tronics,  19,   857-862   (1976). 

[21]     K     P.   D.   Lanyon,  A.   K.   McCurdy,   and  R.  A.   Tuft     "Inf-ared 

3?h°TFeFFP°pnhSe,0f  S!b'CCn  Devices  at  T-15  Microns,"  Record  of 
-^LlEJLPJic^pjta^^  60_63  ('19~HO_or 

[22]     A.   Neugrcscnel     P.   J.   Chen,  S.   C.   Pao,   and  F.  A     Lindholm 

Forward-Bias  Capacitance  and  Current  Measurements  for  Deter 


187 


[24]  M.  A.  Shibib,  F.  A.  Lindholm,  end  F.  Therez,  "Heavily  Doped 
Transparent-Emitter  Regions  in  Junction  Solar  Cells,  Diodes, 
and  Transistors,"  IEEE  Trans.  Electron  Devices,  FD-26,  959- 
965  (1979). 

[25]  E.  L.  Heasell,  "The  Emitter  Efficiency  of  Silicon  Bipolar 

Transistors— An  Unperterted  Band  Model,"  IEEE  Trans.' Electror 
Devices,  ED-26,  919-923  (1979). 

[25]  J.  S.  Blakemore,  Semiconductor  Statistics,  Perg-.mc-n  Press,  New 
Ycrk,  1962. 

[27]  W.  Shockley  and  W.  1".  Read,  "Statistics  of  tne  Recombination  of 
Holes  and  Electrons,"  Phys.  key.,  87,  835  (1952). 

[28]  R.  N.  Hall,  "Electron- Hole  Recombination  in  Germanium,"  Fhvs. 
Rev.,  37,  387  (1952).  ~" " 

[29]  J.  n.  Possum,  F.  A.  Lindholm,  and  M.  A.  Shibib,  ,,The  Importance 
of  Surface  Reccmoinaticn  and  Energy-Bandgap  Narrowing  in  p-n- 
Junction  Silicon  Solar  Ceils,"  IEEE  Trans.  Electron  Devices, 
ED-26,  1294-1298  (1979).       ~~ "  " 

[30]  D.  Redfield,  ''Mechanism  cp  Performance  Limitations  in  Heavily 
Doped  Silicon  Devices,"  Appl .  Phys.  Lett.,  33,  531-535  (1978). 

[31]  H.  U.  Brandhorst,  Jr.,  ''Current  Status  of  Silicon  Solar  Call 
Technology,"  Tech.  Dig.  IEDM,  331-334  (1975). 

[32]     P.   M,   Dunbar  arid  J.   R.   Hauser,   "Theoretical   Effects  of  Surface 
Diffused  Region  Lifetime  Models  on  Silicon  Sclar  Cells," 
Solid-State  Electronics,   20,   697-701    (1977). 

[33]  E.  J.  McGrath  and  D.  H.  Navon,  "Factors  Limiting  Current  Gain 
in  Power  Transistors,"  IEEE  Trans.  Electron  Devices,  ED-24, 
1255-1259  (1977;.        ' 

[34]  P.  Lauwers,  J.  Van  Meerbergen,  P.  Sulteel,  R.  Mertens,  and 
R.  Van  Overstraeter,  "Influence  of  3andgap  Narrowing  on  the 
Performance  of  Silicon  n-p   Solar  Cells,"  Solid-State  Electronics, 
21,747-752(1973).  ~" 

[35]  M.  S.  Adlar,  "Achieveing  Accuracy  in  Transistor  and  Thyristor 
Modeling,"  Invited  Paper,  Tech.  Dig.  IEDM,  550-554 .  1978. 

[36]  J.  Dziewior  and  W.  Schmid,  "Auger  in  Coefficient  for  Highly  Dooed 
arid  Highly  Excited  Silicon,"  Appl .  Pnys.  Lett.,  31,  346-348 
(1977). 


188 


[37]  W.  \L   Sneng,  "The  Effect:  of  Auger  Recombination  on  the  Emitter 

Injection  Efficiency  of  Bipolar  Transistors,"  IEEE  Tra'is.  Electron 
Devices,  ED-22,  25-27  (1975).  "      "    " 

[33]  H.  T.  Weaver  and  R.  D.  Nasby,  "Evidence  for  Bandgap  Narrowing 

Effects  in  Silicon  Solar  Cell  Performance,1'  Tech.  Dig.  IEDH,  693 
1979. 

[39]  C.  T.  Sah  and  C.  T.  Wang,  "Experiments  on  the  Crigin  of  Process 
Induced  Recombination  Centers  in  Silicon,"  J.  AppJ .  Phys. ,  46, 
1757-1776  (1975). 

[40]  J.  Lindmayer  and  J.  F.  Allison,  "The  Violet  Cel7:  An  Improved 
Silicon  Solar  Cell,"  COMSAT  Tech.  Rev.,  3,  1-22  (1973). 

[41]  M.  Klein,  "The  Injection  Efficiency  in  Double  Diffused  Transis- 
tors," Proc.  IRE,  49,  1708  (1961). 

[42]  D.  ?.  Kennedy  and  P.  C.  Murley,  ''Minority  Carrier  Injection 
Characteristics  of  the  Diffused  Emitter  Junction,''  IRE  Trans. 
Electron  Devices,  ED-9:  136-142  (1962). 

[43]  j.  G.  Fossum,  "Computer- Aided  Numerical  Analysis  of  Silicon 
Solar  Cells,"  Solid-State  Electronics,  19,  263-277  (1976). 

T44l  E.  L.  Heasell,  "Injection  Efficiency  in  Bipolar  Devices," 
Tech.  Dig.  IEDM,  310-311  (1978). 

[45]  R.  S.  Popovic,  "On  the  Apparent  Narrowing  of  the  Bar.dgap  in 

the  Base  Region  of  Bipolar  Transistors,"  Solid-State  Electronics, 
22,  248-349  (1979). 

[46]  K.  Seeger,  Semiconductor  Physics,  Soringer-Verlag,  Mew  York, 
1973. 

[47]  G.  L.  Pearson  and  J.  Bardeen,  "Electrical  Properties  of  Pure 
Silicon  and  Silicon  Alloys  Containing  Boron  and  Phosphorous," 

Phys.  Rev.,  75,  365  (1949). 

[48]  N.  A.  Pe-nin,  3.  G.  Zhurkin,  and  B.  A.  Vclkov,  "Influence  cf 
Concentration  of  Donors  and  Acceptors  en  the  Electrical  Con- 
ductivity of  Heavily  Dooed  n-Type  Silicon,"  Sov.  Phys.-- 
Solid  State,  7,  2580  (1965-66). 

[4?]  W.  R.  Thurber,  'Ionization  of  D^oant  Density  in  Silicon," 

p.  13  in  NBS  Special  Publication  400-19.  January  1  co  June  30, 
1975. 

[50]  A.  A.  VoVfson  and  V.  '<.  Subashiev,  "Fundamental  Absorption 

Edge  of  Silicon  Heavily  Doped  with  Donor  or  Accepter  Impurities," 
Sov .  Phys . --Semiconductors ,  1,  327  (1967). 


189 


[51]  6.  A.  Swart?.,  "Low-Temperature  Hall  Coefficient  and  Conductivity 
in  Heavily  Doped  Silicon,"  J.  Phys.  Chem.  Solids,  12,  245  (I960;. 

[52]  P.  W.  Chapman,  0.  N.  Tufte,  J.  D.  Zook  and  D.  Long,  "Electrical 
Properties  of  Heavily  Doped  Silicon,"  .1.  Appl .  Fhys.,  3^,  3291 
(1963).  

[53]  I.  Granacher  and  W.  Czaja,  "Mobility  and  Electron  Spin  Resonance 
Linewidth  in  Phosphorous  Doped  Silicon,"  J.  Phys.  Chem.  Solids, 
28,  231-238  (1967). 

[54]  C.  Yamanouchi,  K.  Mizuguchi,  and  W.  Sasuki,  "Electric  Conduction 
in  Phosphorous  Doped  Silicon  at  Low  Temperatures,"  J.  Phys. 
Soc.  Japan,  22,  859  (1967). 

[55]  F.  Mousty,  P.  Ostoja,  and  L.  Passari,  "Relationship  between 

Resistivity  and  Phosphorous  Concentration  in  Silicon,"  J.  Appl . 
Phys.,  45.  4576  (1974). 

[56]  N.  Kobayashi,  S.  Ikehata,  S.  Kobayashi,  and  W.  Sasaki,  "Specific 
Heat  Study  of  Heavily  P  Doped  Si,"  Sol  id-State  Commun.,  24,  67 
(1977). 

[57]  R.  K.  Sundfors  and  D.  F.  Holcomb,  "Nuclear  Magnetic  Resonance 
Studies  of  the  Metallic  Transitions  in  Dooed  Silicon."  Phys. 
Rev  . ,  136,  A810-820  (1964). 

[58]  M.  N.  Alexander  and  D.  F.  Hclccmb,  "Semiconductor- tc-Metal 

Transition  in  n-type  Group  IV  Semiconductors,"  Rev.  Mod.  Phys., 
40,  815  (1968). 

[59]  C.  Kittel,  Introduction  to  Solid  State  Physics,  4th  ed.,  John 
Wiley  ano  Sons,  Inc.,  New  York,  1971. 

[60]  L.  I.  Rosier,  "Surface  States  and  Surface  Recombination  velocity 
Characteristics  of  S i - S i 0 2  Interfaces,"  IEEE  Trans.  Electron 
Devices,  ED-13,  260-253  (1966). 

[61]  C.  T.  Sah  and  F.  A.  Lindholm,  "Carrier  Generation,  Recombination, 
Trapping,  and  Transport  in  Semiconductors  with  Position- 
Deoendent  Composition,"  IEEE  Trans.  Electron  Devices,  ED-24, 
253-362  (1977).        ""  " 

[62]  P.  M.  Dunbar  and  J.  R.  Hauser,  "A  Study  of  Efficiency  in  low 
Resistivity  Silicon  So"!ar  Cells,"  Solid-State  Electronics,  19. 
95-1 02  ( 1 976) . 

[63]  J.  R.  Hauser  and  P.  M.  Dunbar,  "Performance  Limitations  of 

Silicon  Solar  Cells,"  IEEE  Trans.  Electron  Devices,  ED-24.  305- 
321  (19771. 


190 


[54]  F.  A.  Lindholm,  A.  Neugroshel,  S.  C.  Pac,  J.  £.  Fossum,  and 

C.  T.  Sah,  "Design  Considerations  for  Silicon  HLE  Solar  Cells," 

Record  of  the  13th  Photovoltaic  Specialists  Conf.,  1300-1305 
(1978). 

[65]  A.  Neugroshel,  r.  A.  Lindholm,  S.  C.  Pao,  and  J.  G.  Fossum. 
''Emitter  Current  Suppression  in  a  High-Low-Junction  Emitter 
Solar  Cell  Using  an  Oxide-Charge  Induced  Electron  Accumulation 
Layer,"  Appl .  Phys.  Lett.,  33, "168-170  (1973). 

[66]  P.  A.  lies  and  S.  I.  Soclof,  "Effect  cf  Impurity  Doping  Concen- 
tration on  Solar  Cell  Output,"  Record  of  11th  Photovoltaic 
Specialists  Con*.,  19-24  (1975). 

[67]  M.  Wolf.  "The  Silicon  Solar  Cell  Design  Handbook,"  Record  of 
the  9th  Photovoltaic  Specialists  Conf.,  53-60  (19727!     "" 

[63]  J.  G.  Fossum  and  E.  L.  Burgess,  "High-Efficiercy  p+nn+  Back- 
Sur^ace-Field  Silicon  Solar  Cells,"  Appl .  Phys.  Lett. ,  33, 
23S-24C  (1378). 

[69]  J.  G.  Fossum,  R.  D.  Nasby,  and  E.  L.  Burgess,  "Development  of 
High-Efficiency  p+nn+  Back-Surface-Field  Silicon  Solar  Cells," 
Record  of  the  13th  Photovoltaic  Specialists  Conf.,  1294- 
1299  TWW. 

[70]  S.  S.  Li  and  F.  A.  Lindholm,  "Alternative  Formulation  of 

Generalized  Einstein  Relation  for  Degenerate  Semiconductors," 
Proc.  IEEE,  56,  1256-1257  (1968). 

[71]     I.   V.   Grekhov,   N.   N.    Korobkov,  and  A.   E.  Qtblesk,   "Injection 
Efficiency  of  p-n  Junctions   in  Silicon,"  Sov.   Phys. --Semi con- 
ductors,  9,   4S3-494   (1975). 

[72]  K.  P.  Abdurakhmanov,  Sh.  Mirakhmedov,  and  A.  T.  Teshabaev, 

"Characteristics  of  the  Distribution  of  the  Density  of  States 
in  Heavily  Doped  Silicon,"  Sov.  Phys. --Semiconductors.  12, 
457-459  (1973).  " 

T73]  T.  N.  Morgan,  "Broadening  of  Impunity  Bands  in  Heavily  Doped 
Semiconductors,"  Phys.  Rev.,  139,  A343-343  (1965). 

[74]  H.  K.  Gummel ,  "A  Self-Consistent  Iterative  Scheme  for  One- 
Dimensional  Steady  State  Transistor  Calculations,"  IEEE  Trans. 
Electron  Devices,"  ED- 24,  299-305  (1977). 

[75]  3.  V.  Gokhale,  "Numerical  Solutions  for  a  One-Dimensional 

Silicon  n-p-n  Transistor,"  IEEE  Trans.  Electron  Devices,  ED-17, 
594-602  (1970).  '" 


191 


[75]  5.  D.  Hachtel,  ?,.   C.  Joy, and  J.  I!.  Cooley,  ;'A  New  Efficient 

One-Dimensional  Analysis  Program  for  Junction  Device  Modeling," 
Proc.  IEEE,  60,  36-98  (1972). 

[77]  J.  w.  Slotboom,  "Computer  Aided  Two-Dimensional  Analysis  of 
Bipolar  Transistors,"  IEEE  Trans.  Electron  Devices,  ED-20, 
669-679  (1973).       ~~~ 

[78]  H.  H.  Heimeier,  "A  Two-Dimensional  Numerical  Analysis  of  Silicon 
n-p-n  Transistor,"  IEEE  Trans.  Electron  Devices,  ED-20,  708- 
174  (1974).       "     "'  " 

[79]  0.  Manck  and  M.  L.  Engl,  "Twc-Dimensional  Computer  Simulation 
for  Switching  a  Bipolar  Transistor  Out  of  Saturation,"  IEEE 
Trans.  Electron  Devices,  ED-22.  339-347  (1975). 

[80]  A.  Chryssafis  and  W.  Love,  "A  Computer-Aided  Analysis  of  One- 
Dimensional  Thermal  Transients  in  n-p-n  Power  Transistors," 
Solid-State  Electronics,  22,  249-255  (1979). 

[81]  S.  P.  Gaur,  "Performance  Limitations  of  Silicon  Bipolar  Tran- 
sistors," IEEE  Trans.  Electron  Devices,  ED-26,  415-421  (1979). 

[82]  A.  Phillips,  Jr.,  "Calculations  of  the  Effect  of  Emitter 
Compensation  on  8  and  fT  of  Bipolar  Devices,"  IBM  J.  Res. 
Develop.,  22,  687-689  (1971).  

[83]  H.  J.  DeMan,  "The  Influence  of  Heavy  Doping  on  the  Emitter 
Efficiencv  of  a  Bipolar  Transistor,"  IEEE  Trans.  Electron 
Devices,  ED-18,  833-835  (1971). 

[84]  S.  L.  H.  Wilson,  "The  Emitter  Efficiency  of  Silicon  Transistors," 
Solid-State  Electronics,  20,  71-74  (1977). 

[85]  E.  0.  Kane,  "Thomas-Fermi  Approach  to  Impure  Semiconductor 
3and  Structure,"  Phys.  Rev.,  131,  79-38  (1963). 

[861  P.  A.  Wolff,  "Theory  of  Band  Structure  of  Very  Degenerate 
Semiconductors,"  Phys.. .Rev.,  126,  405-412  (1962). 

[87]  F.  A.  Lindholm  and  R.  W.  Ayers,  "Generalized  Einstein  Relation 
for  Degenerate  Semiconductors,"  Proc.  IEEE,  56,  271-372  (1968). 

[88]  0.  Madelung,  Handbuch  der  Physik,  20,  Springer,  Berlin,  1957. 

[89]  J.  McQougall  and  E.  C.  Stoner,  "The  Computation  of  Fermi-Dirac 
Functions,''  Phil.  Trans.,  A  237,  67-104  (1938). 

rgol  W.  Ehrenberg,  "The  Electric  Conductivity  of  Simple  Semicon- 
ductor," Proc.  Phys.  Soc,  A  63,  75-76  (1950). 


191 


[91]  ?.  T.  Landsberg,  R.  W.  MacKay  and  A.  D.  McRonalri,  "The  Parameters 
of  Simple  Excess- Semiconductors,"  Proc.  Phys.  Soc,  A   476-480 
(1950). 

[92]  R.  B.  Dingle,  "The  Fermi-Dirac  Integrals, "  Appl .  Sci.  Res., 
36,  225-239  (1955-57). 

T93]  R.  E.  Eattocletti,  "Polynomial  Approximation  of  Fermi  Integrals," 
Proc.  IEEE,  53,  2162-2153  (1965). 

[94]  E.  L.  Jones,  "Rational  Chebyshev  Approximation  of  the  Fermi- 
Dirac  Integrals,"  Proc.  IEEE,  4,  70S-7C9  (1966). 

[95]  A.   B.  Joyce  and  R.  W.  Dixon,  "Analytic  Approximation  for  the 
Fermi  Enerqy  of  an  Ideal  Fermi  Gas,"  Appl .  Phys.  Lett- ,  31, 
354-356  (1977). 

[96]  K.   G.  Nilsson,  "Empirical  Approximations  for  the  Fermi  Energy 

in  a  Semiconductor  with  Parabolic  Bands,"  Appl.  Phys.  Lett..  33, 
653-654  (1973). 


BIOGRAPHICAL  SKETCH 

Muhammed  Aymar,  Shi  bib  was  born  in  Damascus,  Syria,  on  February 
14,  1953.  He  completed  an  advanced  high  school  degree  in  mathematics 
and  physics  at  the  International  College,  Beirut,  Lebanon,  and 
graduated  in  June  1972  with  high  distinction.  In  June  1975,  he  received 
the  B.S.  degree  from  the  American  University  of  Beirut,  Eeirut,  Lebanon. 
His  major  field  was  in  physics,  and  he  was  or.  the  American  University 
of  Beirut  Dean's  Honor  List,  He  started  graduate  studies  in  electrical 
engineering  at  the  University  of  Florida  in  September  1975,  dni   he 
graduated  in  December  1976  wi bh  an  M.S.  and  in  December  1979  with  a  Ph.D. 

His  master's  researcn  involved  computer-aided  analysis  of 
Integrated  Injection  Logic  (I^L).  Since  1975.  he  has  been  a  graduate 
research  assistant  and  a  graduate  research  associate. 

His  research  interests  are  in  solar  calls,  semiconductor  and 
solid-state  physics,  and  semiconductor  device  modeling. 

He  is  a  member  of  IEEE,  the  American  Physical  Society,  and 
the  Society  of  Physics  Students. 


193 


I  certify  that  I  have  read  this  stud>  and  that  in  my  opinion  it 

conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 

adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


Fredrik  A.  Lindholm,  Chairman 
Professor  of  Electrical  Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  oninion  it 

conforms  to  acceptable  standards  of  scnolarly  presentation  and  is  fully 

adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


Jerry/  G.  Fossum,  Co-Chairman 
Associate  Professor  of  Electrical 
Eiiginearing 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


Arnost  Neugroscnel 
Associate  Professor  of  Electrical 
Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 

conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 

adequate,  in  sccp2  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


SWehg   S.  Li  S 

Professor  of  Electrical  Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 

conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 

adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


V 


77" 


Chanes'F.  Hooper,  Jr, 
Professor  of  Physics 


This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the  College 
cf  Engineering  and  to  the  Graduate  Council,  and  was  accepted  as  partial 
fulfillment  of  the  requirements  for  the  degree  of  Doctor  of  Philosophy. 


December  1979 


--^. 


I/*-/ 


'% 


■^O^L-^^i 


an,  Co'i^ge  of  Engineering 


.<Z2=^. 


Dean,  Graduate  School 


UNIVERSITY  OF  FLORIDA  . 

■H 


